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ABSTRACT 

We present an analysis of the X-ray point source populations in 182 Chandra 
images of galaxy clusters at z > 0.1 with exposure time >10 ksec, as well as 44 
non-cluster fields. Analysis of the number and flux of these sources, using a detailed 
pipeline to predict the distribution of non-cluster sources in each field, reveals an 
excess of X-ray point sources associated with the galaxy clusters. A sample of 148 
galaxy clusters at 0.1 < z < 0.9, with no other nearby clusters, show an excess of 230 
cluster sources in total, an average of ~1.5 sources per cluster. The lack of optical 
data for these clusters limits the physical interpretation of this result, as we cannot 
calculate the fraction of cluster galaxies hosting X-ray sources. However, the fluxes 
of the excess sources indicate that over half of them are very likely to be AGN, and 
the radial distribution shows that they are quite evenly distributed over the central 
1 Mpc of the cluster, with almost no sources found beyond this radius. We also use 
this pipeline to successfully reproduce the results of previous studies, particularly the 
higher density of sources in the central 0.5 Mpc of a few cluster fields, but show that 
these conclusions are not generally valid for this larger sample of clusters. We conclude 
that some of these differences may be due to the sample properties, such as the size 
and redshift of the clusters studied, or a lack of publications for cluster fields with 
no excess sources. This paper also presents the basic X-ray properties of the galaxy 
clusters, and in subsequent papers in this series the dependence of the AGN population 
on these cluster properties will be evaluated. 

In addition the properties of over 9500 X-ray point sources in the fields 
of galaxy clusters are tabulated in a separate catalogue available online or at 
www.sc.eso.org/ ^rgilmour. 

Key words: galaxies: active , galaxies: clusters: general , X-rays: galaxies , X-rays: 
galaxies: clusters 



1 INTRODUCTION 

Studies of AGN host galaxies are key to understanding the 
physical mechanisms which trigger AGN activity, and govern 
the fueling rate of the central black hole. An important part 
of these studies is the external environment of the galaxies, 
which has long been kno wn to have a significant l ink with 
the galaxy properties fe.g lHubble fc Hu masonlll93lh. Corre- 
lations such as the morpholog y-density (Drossier 198^) and 
star- formation - density (e.g. Gome z et al.ll2003i ') relations 
are evidence of the significant transformations that are as- 
sociated with galaxy clusters. If AGN activity is influenced 
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by, for example, galaxy mergers or gravitational disruption 
of the host galaxy, then the number of AGN would also differ 
between galaxy clusters and the field, and within the cluster 
itself. 

The first evidence of a link between AGN activity and 
environment came in 1978, when iGisleil (|l978l ) found a lack 
of emission-line galaxies in galaxy clusters relative to the 
field, which was confirmed by Dressier, Thompson & Sheet- 
man (1984) . More recently, large optical surveys have been 
used to identify the properties of a significant number of 
host galaxies, and hence compare AGN activity in the same 
type o f host galaxy but different environment. iMiller et al.l 
l|2003l ) find t hat optical AGN act ivity is independent of envi- 
ronment, but I Wake et al.l l|2004l ) show that the level of clus- 
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taring depends on AGN luminosity. iKauffmann et alj (|2003l ) 
find that AGN wit h strong [Q UI] emission avoid areas of 
high galaxy density. iBesd l|2004 ) investigate the radio prop- 
erties of these AGN, and find that the fraction of galax- 
ies with radio-loud AGN increases dramatically with local 
galaxy density, but that all of these AGN have low [OIII] 
emission and so may not be seen in optical surveys. 

Arguably the least biased method of detecting AGN 
currently is to use X-ray images, which have the added ad- 
vantage that the vast majority of point sources are AGN. 
Not long after optical surveys identified a lack of emission- 
line galaxies in clusters. X-ray surveys began to find a 
surprisingly hig h number of point sources in fields with 
galaxy clusters (iBechtold et all Il983l : iHenrv fc Br"ieil Il99ll : 
iLazzati et al.l 1 19981 ) . With the advent of the Chandra X- 
ray telescope, with sub-arcsecond point sources, such stud- 
ies were repeated for other clusters, with a range of results. 
Significant overdensities of point sources have been found 
in a n umber of fields wi th gala xy clusters at m oderate red- 
shifts llCappi et al.ll200ll . z=0. 5;lMartini et al.| [2b02. z=0.15; 

"2001 z=0.32; Ijohnson et al.l |2003|, z= 0.83; 

< z < 0.31; iD'Elia et al.ll2004 



Molnar ot al 

Martini ct al. 2006. 0.05 ^ 

z=0.5) and gr oups (|Jeltema et al.ll200ll . 0.2 < z < 0.6), but 
iMolnar et al.l also found a z=0.5 cluster without a signifi- 
cant overdensity. More re cently, studies of signifi cant sam- 
ples of galaxy cluster s bv ICappelluti et al.l ( 2005 ) (10 clus- 
ters, 0.24 < z < 1.2). lRuderman fc EbeUngI boOSi) (51 clus- 
ters, 0.3 < z < 0.7), and Branches) et al. 1 20071 ) (18 clusters, 
0.25 < z < 1.01) have all found significant overdensities of 
point sources over the full sample, but not neces sarily in all 
indiv idual fields. However, the ChaMP project (|Kim et al.l 
I2OO4I ) find no difference in the number density of sources 
in fields with z > 0.3 clusters compared to those without 
clusters. 

The calculated number of AGN per cluster varies signif- 
icantly in these samples, even taking into account the differ- 
ent depths of the observations and the statistical variance 
in the number of background sources in each image. This 
is to be expected as the number of AGN per cluster is, of 
course, related to the cluster properties, such as the number 
of possible AGN host galaxies. However, due to the lack of 
optical data it is hard to draw any conclusions from these 
samples as to how, if at all, the cluster environment affects 
the AGN population. Optical imaging and spectroscopy can 
identify the X-ray detected AGN in the cluster, rather than 
relying on statistical background subtraction, and also re- 
veal the distribution and number of normal cluster galaxies. 
This would allow the calculation of the fraction of cluster 
galaxies which host AGN, as a function of cluster radius 
or cluster size for example, shedding light on the physical 
mechanisms affecting AGN. On the other hand obtaining 
optical data for a large sample of galaxy clusters is time 
consuming, and any strong trend should be visible in the 
distribution of X-ray point sources in a large sample of clus- 
ter fields. A rough estimate of the cluster galaxy population 
can also be made from the extended X-ray gas. An analysis 
based purely on X-ray data can therefore be very useful, but 
is clearly inferior to a full spectroscopic analysis of a large 
sample of clusters, with X-ray and optical data. 

Such a study was started bv iMartini et al.l l|2002l ). us- 
ing optical data to identify X-ray detected AGN in galaxy 
clusters. Martini, Mulchaey & Kelson (2007) have confirmed 



optically that the fraction of galaxies hosting X-ray de- 
tected AGN does differ significantly between the eight clus- 
ters in their sample, implying that the cluster properties 
affect the number of AGN. Their spectroscopic data con- 
firms between 2 and 10 AGN per cluster, with a range in 
AGN fractions that cannot be explained by Poissonian vari- 
ations. The mean fraction for the whole sample is 5 per cent 
of galaxies with Mr < -20 hosting AGN with Lx > W*^ 
erg sec~^. The wide variation in the number or fraction of 
X-ray AGN is also found when other stu dies with spectro- 
scopic data are compared. For example iFinoguenov et ahl 
(2004) find only one confirmed X-ray AGN, with luminos- 
ity ~ lO*'^ erg sec"'^, in a 1.8 square d egree survey of the 
centre of the Coma cluster (z=0.02), but lDavis et al.l (|2003l ) 
find between 3 and 5 AGN in a cluster at z=0.08, giving a 
fraction o f 4 p er cent in agreement with the mean value of 
Mar tini et al.l 

iMartini et al.l also find tentative evidence that the clus- 
ters with higher AGN fraction have lower redshift, lower 
velocity dispersion, higher substru cture and lower Butcher- 
Oemler (|Butcher fc Oemleij [l984l ) fraction, but due to the 
small size of the sample (8 clusters) it is not clear from 
this sample which, if any, of these factors is affecting the 
AGN activity. A higher AGN fraction at high redshift 
would be expected from the fie ld evolution of AGN, but 
no strong evolution i s found by iBranchesi et al. I (l2007l) or 
iRuderman fc EbeUnd (|2005l ). In contrast. ICappelluti et al] 
I 2005h find some evidence for an increase in AGN with red- 
shift, and lEastman et al.l (|2007l ) conclude that the increase 
of bright AGN in clusters is up to 20 times greater than in 
the field between 2 ~ 0.2 and z ~ 0.6. 

If AGN are triggered by galaxy mergers then clus- 
ters with lower velocity dispersions would be expected to 
have hig her AGN fr actions, as they have a higher merger 
rate. Po pesso fc Bivi ano (200f|) show that this is indeed the 
case for optically detected AGN, with clusters with high 
velocity dispe r sions having lower AGN fractions. However 
IMartini et all l|2007l ) find that the velocity distribution of 
AGN in eight clusters is not significantly different from the 
velocity distribution of the non-active cluster galaxies. 

Galaxy mergers are also more common in the outskirts 
of clusters, so the projected radial distribution of AGN 
should be diffe rent from the host galaxies if mergers cause 
AGN activity. IMartini et al.l find no evidence for AGN to 
lie in galaxies in the outskirts of the cluster compared to 
bright cluster galaxies, in fact on the contrary, they find ev- 
idence that galaxies with luminous AGN (> 10'''^ erg sec~^) 
are more centra lly clustered than the gen eral population. 
IBranchesi et al.l and iRuderman fc Ebelind also present ev- 
idence that most AGN ar e found within the ce ntral 0.5 
Mpc of the clusters, which IRuderman fc Ebelind attribute 
to tidal encou nters with the central galaxy. In contrast, 
Ijohnson et al.l (,2003) find the excess AGN in a cluster at 
z=0.83 lie between 1 and 2 Mpc from the cluster centre. 
iGilmour et ai] |20o3) also find that AGN avoid the densest 
areas of a supercluster at z=0.17. 

The wide range of results from the current studies may 
partly be due to the large number of variables which can 
affect the AGN fraction in clusters, and partly due to sta- 
tistical fluctuations in the small number of AGN found 
per cluster. In addition, studies which count X-ray point 
sources, without optical confirmation of cluster member- 
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ship, are limited by field-to- field variation in the clustering 
of background sources (e.g. iGiUi et"aLll2005l ). Careful data 
reduction is also required, as detecting AGN against the 
extended intra-cluster medium and accounting for the vari- 
ations in the point spread function (PSF) is important in 
determining the expected number and distribution of back- 
ground AGN in the field. Different treatments of these vari- 
ables may account for the differences between the results of 
the ChaMP project and the surveys which deliberately tar- 
get AGN in cluster emission. In order to understand how 
many AGN are in clusters with different properties, and to 
get significant statistics, a large sample of galaxy clusters 
and non-cluster fields are required. Even in a such a sample 
the conclusions are limited by the lack of optical data which 
means that individual cluster AGN cannot be identified. We 
may, however, expect to see any strong trends in the aver- 
age number of AGN in clusters of a given type, or at a given 
epoch. As shown later in this paper, at least five fields are 
required in each subsample to remove the effects of cosmic 
variance in the background distribution. Larger samples are 
likely required in order to get statistically significant results. 



2 OUTLINE AND METHOD 

The large number of observations of galaxy clusters in the 
Chandra archive provide an excellent basis for investigating 
the prevalence of AGN in galaxy clusters. By comparing the 
point source distribution in 'blank field' observations with 
that found in cluster observations, the number, flux and ra- 
dial distribution of the sources associated with the cluster 
can be determined statistically. In addition the Chandra field 
of view allows AGN to be detected accurately up to 8 ar- 
cmin from the centre of the field. This method has been used 
in the past to investigate small samples of galaxy clusters, 
but these contain significant errors due to field-to-field vari- 
ance. This study is a significant advance over previous stud- 
ies, both in size and methodology. By analysing 182 galaxy 
clusters the statistical variance seen in the smaller studies is 
significantly reduced, and the properties of the cluster AGN 
population can be identified. Furthermore a sample of this 
size can be split into sub-samples and still produce signif- 
icant results. The dependence of the AGN population on 
cluster redshift, mass (estimated from the X-ray luminos- 
ity) and morphology can therefore be found. This analysis 
requires careful data reduction and modelling of the sensitiv- 
ity of each observation to point sources, which varies across 
the image. This was performed using an automated pipeline 
developed for this purpose. 

The key steps in investigating the point sources in the cluster 
observations are as follows: 

• Observations of galaxy clusters with published redshifts 
> 0.1 and 'blank' fields from the Chandra archive are se- 
lected and reduced. 

• Each image is visually inspected to insure that the clus- 
ter is detected at the expected location, and that the image 
does not contain multiple clusters. The luminosity of each 
cluster is found and used to estimate the effects of gravita- 
tional lensing on the background sources. The cluster lumi- 
nosity and assigned morphological class (see Section I3.2.ip 
also allow a later comparison of the AGN content of clusters 
as a function of cluster properties. 



• Point sources are identified in the fields, and their prop- 
erties are calculated. 

• For each observation a 'flux-limit map' is produced, 
showing the detection sensitivity at each point on the im- 
age. This accounts for the detector response, size of the PSF, 
and the level of background emission, particularly from the 
intra-cluster medium. 

• The Log N(> S) - Log S distribution (where N is the 
number of sources and S is the flux) is calculated for each 
blank and cluster fleld, taking into account the sky area 
sensitive to sources of each flux value. 

• The radial distribution of sources, as a function of dis- 
tance from the cluster centre, is calculated. A predicted 
radial distribution, assuming no cluster AGN, is produced 
from the blank field source distribution and the flux-limit 
map. 

• The effects of gravitational lensing of background X-ray 
sources by the galaxy cluster are modelled, and the Log N(> 
S) - Log S distributions and predicted radial distributions 
are corrected for this effect. 

The number, flux and radial distribution of the X-ray 
sources in clusters can then be determined statistically by 
comparing the actual results for each cluster with the pre- 
diction, which assumes that no cluster AGN exist. Section[3] 
describes the data reduction and sample selection, for both 
cluster and blank fields. Section |4] explains the source de- 
tection, and Section [5] the model for producing a predicted 
distribution. Section|5]explains the first results of this study. 
Further results will be published in an accompanying paper. 



3 INITIAL DATA REDUCTION AND SAMPLE 
SELECTION 

There were around 700 imaging observations marked as 
'Clusters of Galaxies' in the Chandra archive in mid-2007. 
However the majority of these are not valid for this study, 
for a range of reasons. In order to determine which observa- 
tions are useful it is necessary to first reduce the data, as only 
then can the reality and position of the cluster be found. An 
initial sub-sample of these observations was therefore put 
through the first stage of the automated pipeline before the 
final sample was defined. This sub-sample contained all ob- 
servations with published redshift > 0.1 and exposure time 
> 10 ksec. The details and reasons for these cuts, and the 
further restrictions applied to produce the final sample, are 
described in Section [3.21 

3.1 Data reduction 

To reduce the initial cluster and blank field samples, an auto- 
mated pipeline was developed, using a range of CIAO tools 
and other programs. This ensured that the reduction was 
uniform, and allowed the whole sample to be reduced effi- 
ciently. In order to obtain the maximum number of sources 
around each cluster, all four ACIS-I chips were used for ob- 
servations focused on the ACIS-I array, and the three chips 
nearest to the aim-point were used for ACIS-S observations 
(or less if not all were turned on). Due to the off-axis degra- 
dation of the PSF the other chips were not investigated as 
the errors become too large. Parts of the selected chips were 
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later excluded as the analysis was restricted to a maximum 
radius from the aim-point (see Section 15. ip . 

For each observation the data were re-reduced from the 
level 1 event list using standard CIAO 3.0.1 tools. The 
fix-hatc^ script was used to check and correct the astrom- 
etry for systematic aspect offsets. A time dependent charge 
transfer inefhciency correction was applied in all observa- 
tions taken after 29 January 200(13. A new level 2 file was 
created, using CALDB 2.26 to correct for the degradation of 
the QE. The data were filtered for standard grades, status=0 
and the default good time interval (GTI). Further GTI fil- 
tering was performed for each chip by manually masking 
the brightest sources and filtering for count rates more than 
3-sigma above the quiescent value. CCD 8 was destreaked 
using the standard tools, and the data were filtered for bad 
pixels. Finally, the data were filtered for energies between 
0.5 and 8 keV to allow better detection of AGN. 

3.2 The cluster sample and cluster properties 

An initial sample of observations from the 'Clusters of 
Galaxies' category, with (probable) redshift > 0.1 and ex- 
posure time > 10 ksec was selected for this project. Lower 
redshift clusters were excluded as only the central regions 
would be covered by the Chandra image. For example a 
z = 0.1 cluster observed with the ACIS-S array would be 
observed to a radius of at least 220kpc in all directions, and 
with ACIS-I this increases to at least 440kpc. The maximum 
radius covered is significantly larger than this, as the cluster 
is rarely placed in the centre of the array. 

Regardless of the cluster selection criteria the sample 
will be heavily biased, as clusters are selected depending on 
the requirements of the observer. In particular the sample 
will be biased towards relaxed clusters , which are used t o 
constrain cosmological parameters fe.g. lAllen et al.ll2004al ). 
and rich, highly disturbed clusters, used to study cluster 
mergers. Other observations were searches for cluster emis- 
sion. By examining the proposal abstracts clusters were ex- 
cluded if they were deliberately targeted due to their lensing 
of background QSOs. The remaining biases in the sample se- 
lection were parameterised as far as possible by examining 
the X-ray properties of the clusters, and taken into account 
later in the analysis. 

The cluster redshift was determined from sources in the 
NASA Extragalactic Database (NED). Observations were 
selected which have a confirmed galaxy cluster, cD galaxy, 
QSO or galaxy overdensity at z > 0.1 within 5' of the aim- 
point. The archive was examined up to April 2007 and 192 
targets were selected, of which 34 were observed on more 
than one occasion (with the same detector array). A small 
number of cluster observations fulfilled the above criteria 
but were not suitable for the pipeline due to non-standard 
settings which were not easily incorporated into the data 
reduction. 

Five properties were evaluated for each cluster field - 
the reality, number of clusters, centre, morphology and lu- 
minosity. The first two were used to reject clusters with no 

^ see |http: / / cxc.harvard.edu/cal/ASPECT/fi x_offset / fix_offset.cgi 
^ after the ACIS focal plane temperature was 
lowered |http://cxc.harvard.edu/cal/Acis/CaLprods/tgain/index.html| 



X-ray emission, which are possibly not true clusters, and 
fields with multiple clusters at different redshifts. The X-ray 
position of the cluster is important as the AGN distribution 
may depend on cluster radius, and many optically discov- 
ered clusters have poorly defined centres. The latter two 
properties are evaluated in order to determine whether the 
cluster properties affect the number or distribution of AGN. 
The luminosity also provides an estimate of the mass, which 
can be used to correct the predicted source counts for each 
image for the effect of gravitational lensing, as described in 
Section 15.31 



3.2.1 Cluster reality and spatial properties 

The morphology, centre and reality of each cluster, and the 
number of clusters in each field were determined by exam- 
ining the 0.5-8 keV images and the smoothed background 
images (with point sources removed, see Section [5.1|l . The 
cluster centre was taken to be the peak of the smoothed 
background image. In the few cases where the cluster con- 
sisted of two peaks of similar brightness, the mid-point was 
chosen. 

The reality and morphology were determined by eye, us- 
ing the full and smoothed images (described in Section 15.11) . 
The morphological classifications clearly involve a certain 
amount of subjectivity, but they will be sufficient to iden- 
tify the most disturbed clusters. The following categories 
were used, and are illustrated in Figure [1] 

0. No cluster emission visible against the background 
fiuctuations. 

1. One relaxed cluster. It may be elliptical or have edge 
structure, but not enough to fall into another category. 

2. One disturbed cluster. The disturbance must be such 
that the cluster is clearly not simply elliptical or an asym- 
metric ellipse, and must be joined to the cluster by visible 
emission. 

3. Merging cluster. A double-peaked system, with a sub- 
cluster with peak emission (in the smoothed image) > 20 per 
cent of the main cluster peak, joined to the main cluster by 
visible emission or at the same redshift. 

4. Two clusters. A second cluster with peak emission 
> 20 per cent of the main cluster peak but not clearly asso- 
ciated with it. 

Ic, 2c, 3c. As 1, 2 and 3, but with a small contami- 
nating cluster or group in the field of view. The secondary 
emission must have a peak value of < 20 per cent of the main 
cluster peak and not be clearly associated with it. The con- 
tamination can also be an optically confirmed cluster with 
no X-ray emission, which is in or very near the field of view, 
as described below. 

The images were assigned a category by two observers, 
who were broadly in agreement as seen in Table [T] The few 
discrepancies are mainly due to small contaminating clusters 
which may be background fluctuations or faint undetected 
point sources, and cases of uncertainty over the degree of 
disturbance. The morphologies of Observer 1 (the first au- 
thor) were adopted as they are slightly more conservative. 

The optical data (from the NED) for the cluster fields 
were used to check for optically detected clusters in or near 
to the field of view (within 15'from the cluster centre). Three 
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Morphology Ic Morphology 2 




Morphology 3 Morphology 4 



Figure 1. Examples of cluster morphology categories. The im- 
ages are around 18 arcmin x 12 arcmin. 



this is likely to be minimal in most cases). The final sample 
therefore consists of the 148 observations with morphology 
class 1, 2 or 3 and 0.1 < 2 < 0.9. 

The twenty weakly contaminated cluster fields (Ic, 2c 
and 3c) were included in a second sample, as the fields may 
still be of interest. The ten fields which clearly contained a 
second cluster (type 4) were rejected from the rest of the 
analysis. 

The fourteen 2 > 0.9 cluster observations with secure 
redshifts were placed in a third sample, regardless of the re- 
ality or extent of their emission, as at this redshift range 
almost all of the observed clusters are centred on active 
galaxies, and often the extended emission may be too faint 
to detect or be contaminated by AGN jets Some of these 
objects are better classified as proto-clusters, so they are 
analysed separately from the rest of the sample. 

The final cluster fields, split into the above categories, 
are described in Tabled 



Observer 1 
Ic 2 2c 



3c 



O 





1 

Ic 
2 

2c 
3 
3c 
4 



10 



115 4 2 

13 1 
4 14 

1 1 



13 



1 

10 



Table 1. Morphology classes assigned to 0.1 < z < 0.9 clusters 
by two observers (using X-ray and optical data). 



fields were moved from morphology class 1 to Ic as they con- 
tained optically confirmed clusters at a significantly difi'erent 
redshifts from the main cluster. The images were compared 
to the NED to check that the detected peak in the X-ray 
emission corresponded to the location of the galaxy cluster 
as given in the NED; clusters were accepted if the centre of 
the X-ray emission was within 3' of the NED object. Two 
cluster observations were removed as their optical position 
was > 3' from the observed X-ray peak, and therefore the 
optical redshift may not apply to the X-ray detected cluster. 
For clusters that had more than one redshift measurement 
in the NED the cluster redshift was accepted if 5z/z < 0.1. 
Otherwise the literature was examined in detail to determine 
the most accurate redshift - these clusters are flagged in Ta- 
ble [2] In addition one bright cluster has a slightly revised 
redshift due to an iron emission line in the X-ray spectrum, 
as described in Section [3.2.31 



3.2.2 The final cluster sample 

The 192 cluster fields were split into samples depending on 
their properties. The final cluster sample contains only un- 
contaminated confirmed clusters to ensure that the analysis 
is not affected by additional clusters in or near the field of 
view, which could also contain AGN and may contribute to 
the lensing of background AGN (see Section 15.31 although 



3.2.3 Cluster luminosities and temperatures 

To compare clusters at different redshifts, luminosities need 
to be found in the same rest-frame band for each cluster. A 
spectrum was extracted from the level 2 data for each cluster 
and fit with a thermal model, which was then evaluated in 
the given band. The following analysis was applied to all 
cluster observations, but is only truly valid for clusters with 
morphology classes 1-3 as listed in Table [2] 

Spectra were extracted from the 0.5-8 keV band data 
(to simplify the data reduction) from circular apertures cen- 
tred on the cluster centre. The chosen aperture included 
~ 99.5 per cent of the cluster counts, and the background 
spectrum was taken from an annulus with radii of 1.1 and 
1.49 times the cluster radius. Point source regions were 
subtracted and the regions containing the brightest point 
sources were enlarged if necessary, to ensure that they did 
not contaminate the cluster emission. Areas of bad or no 
exposure were also removed. Response functions were cal- 
culated for the central region of the cluster aperture, rather 
than finding a weighted response over the full aperture, due 
to the time required for the latter. Tests on three clusters 
found that the difference in fiux for a single central response 
file compared to that for the full region was < 2 per cent, 
which is negligible compared to the errors in the model. 

Spectra were fitted using an a bsorbed Raymond- 
Smith model (|Ravmond fc Smithll 19771 ) in XSPEC vll.3.1, 
binned to a minimum of 25 counts per energy interval. The 
gala ctic neutral hydrogen den sity was fixed at the local val- 
ues l|Dickev fc Lockman|[l990l ). and the redshift fixed to the 
value in Table [2] For clusters with more than one observa- 
tion the multiple spectra were fit simultaneously. The model 
errors are underestimated as they do not take into account 
errors due to taking the calibration of the central pixel only. 
To get a better measure of the accuracy of the luminosities, 
the spectra of clusters which were observed twice were fit- 
ted individually, and the difference between the luminosities 
were found to be less than 0.1 dex at all fluxes. The observed 
luminosities of the clusters using this method are given in 
Table [2] These ge nerally match the C handra luminosities in 
the literature (e.g. lEbeling et aLll2007f ) to within 10 per cent. 
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3.3 Blank fields 

It is necessary to have a control sample of blank fields in or- 
der to calculate the expected distribution of point sources in 
each cluster observation due to foreground and background 
objects (which will be referred to as 'background sources', 
although they may in fact be in the foreground). To avoid 
biases due to large scale structure and statistical variance 
due to low counts it is desirable to have as large a sample 
as possible of blank fields. 

Many X-ray surveys of 'blank' fields have been con- 
ducted in order to study the general X-ray source popu- 
lation. Some of these observations were selected from the 
archive, and reduced with the pipeline to ensure consistent 
data reduction. Fields that contained galaxy clusters which 
were discovered independently of the blank field observation 
were removed, but fields with serendipitous cluster detec- 
tions were retained. Individual pointings were selected so as 
to maximise the sky area and match the depth of the cluster 
observations. 

In addition to the 'true' blank fields, observations of 
high-redshift (z > 2) quasars or radio galaxies were also 
used. These were added to increase the sample size, and 
hence reduce the errors due to low source counts (particu- 
larly at high fluxes). In addition, all of the blank field ob- 
servations used the ACIS-I detector, so observations using 
ACIS-S were required to test for differences due to the de- 
tector. The flelds all have observation times > lOksec and 
redshifts in the NED. In most of these flelds the QSO is vis- 
ible at the aim-point and it is possible that there are extra 
sources at the redshift of the QSO due to either clustering 
or lensing (these will be rare as the observations are shallow 
and the target QSOs very distant). In all images a circle 
of radius 25" was removed from around the aim point, as 
this radius excludes all other objects identified in NED at 
the QSO redshifts (with the exception of one field, which 
was rejected). These regions were excluded from the analy- 
sis using the masks described in Section [5. II One field had a 
further region removed due to a rare serendipitous detection 
of a nearby galaxy with resolved point sources. 

Once the data were reduced, the blank field number 
counts were checked to ensure that including the high red- 
shift QSO fields does not bias the background (as explained 
in the Appendix). The final sample of blank fields consists 
of 22 true blank fields and 22 QSO flelds, which are listed 
in Table 



4 POINT SOURCE DETECTION AND 
PROPERTIES 

4.1 Source detection 

Images were made using unbinned data and exposure maps 
(in s~^cm^^) were made assuming the sources have a photon 
index of F = 1.7, typical of unobscured AGN at the sample 
flux limits (see for example flgure 3 of iTozzi et al.l . 2001). 
Tests on a few images showed that changing the spectral 
index to other realistic values does not significantly change 
the sources detected or their significa nces. Sources were de - 
tected using the wavdetect package jFreeman et al.|[2"002l l . 
with wavelet scales of 1,2,4,8 and 16 pixels and a signifi- 
cance threshold of 10"''. Tests using different wavelet scales 



suggest that ^ 1 per cent of sources are close enough to 
be missed by using these scales, but would be detected us- 
ing scales separated by y/2. The source list output from 
WAVDETECT was examined by eye to remove detections of 
the extended cluster emission. A montecarlo simulation of 
the source detection (Appendix A.l) shows that very few 
sources are missed by wavdetect, even accounting for the 
rapidly varying background in regions near the cluster cen- 
tres. 

Many clusters and blank flelds were observed more than 
once, and where possible in these cases data from up to 
three observations were merged before the sources were de- 
tected, to give far deeper images and maximise the number 
of sources. Images were only merged if the same detector 
(ACIS-I or ACIS-S) was used. The process is similar to that 
for single observations with the following additions:- 

The astrometry of the images was adjusted using 
the align.evt routin^fl as, even after correcting the aspect 
flies, small offsets often exist between images. Images were 
matched using sources detected in the central 4 arcmin of 
each image, where the PSF is smallest. Individual images 
were made, and exposure maps created for each observa- 
tion. A combined image and combined exposure map were 
computed. 

WAVDETECT determines whether a source is real based 
on the source extent and the size of the PSF, which is com- 
plex for merged images with different aim points. In this case 
the combined PSF size at each point was calculated by com- 
bining the PSF sizes of the individual images, weighted by 
their exposure map, and calculating the 3-sigma encircled 
energy size of the resulting source. These PSF sizes were 
input into WRECOrjf] in order to give detections and sizes 
that are comparable to the standard wavdetect results for 
single images. 

As an illustration of this technique. Figure [2] shows the 
sources detected in a combined image of MACS J1149+22. 
The combined exposure map and expected PSF distribution 
are also shown. 

4.2 Source properties 

It is important to determine the source properties very ac- 
curately, as only a few percent of the detected sources are 
likely to be cluster sources. The observations have a wide 
range of exposure times, and the source sizes also change 
with off-axis angle, and small errors in the determination 
of the source properties could therefore wipe out any sig- 
nal from the cluster sources, or introduce biases with, for 
example, redshift or exposure time. 

Because of the need for high accuracy the reality and 
properties of the wavdetect detected sources were re- 
determined using more stringent criteria. This is also neces- 
sary in order to create an accurate model of detection prob- 
ability, taking into account the effect of the variable back- 
ground in cluster flelds, as described in Section [S] wavde- 
tect outputs were used to determine source positions and 
sizes, but other properties of the sources, such as counts and 

3 ALIGN_EVT vl.6, written by Tom Aldcroft 

* Using a more flexible version, kindly provided by Peter Freeman 

(private communication). 
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Figure 2. Source detection inputs and outputs for MACS 
J1149+22. The top panel shows the exposure map (left, higher 
exposure is white) and calculated expected PSF size distribu- 
tion (right, smaller PSF is whiter) for this observation. The lower 
panel shows the sources detected on the combined image, using 
these inputs, and a enlarged portion of the image, covering an 
area where the expected PSF size varies rapidly. The full image 
is ~24 arcmin across, and the pixel size is ~0.5 arcsec. It is clear 
that the input PSF size distribution, combined with the detection 
power of WAVDETECT, accurately finds the true source centre and 
extent. Some detected sources could be background fluctuations, 
but these are later removed as described in Section |4] 



significance, were re-calculated. The positions and proper- 
ties of some of the significant sources detected in one field 
are listed in Table H) which also contains the web address of 
the full source list for all cluster fields. 

In order to maximise the signal-to-noise for individual 
sources, the wavdetect source sizes were used to deter- 
mine a circular aperture size for each source (as the current 
Chandra PSF models are only measured at a few radii). The 
aperture had radius 1.2 x rmax, where rmax is the semi-major 
axis of the wavdetect 3-sigma ellipse. Extensive testing 
showed that this radius of aperture maximised the signal- 
to-noise for the sample whilst minimising the missed source 
counts. 

In deep X-ray images many sources either overlap with 
other sources or with areas of bad exposure such as chip 
gaps. Pixels in the source aperture were rejected if they were 
within the aperture of another source, or had exposure be- 
low 10 percent of the median value for the source. These 
pixels were replaced by their refiection on the opposite side 
of the aperture if possible, or otherwise with other pixels 
from the same radius. Around 3 per cent of the sources de- 
tected required some degree of correction, and for < 0.5 per 
cent of sources the correction is only accurate to within a 
factor of ~2 due to the large correction area. 

For each source the mean background count rate per 
pixel was calculated in an annulus of area 10000 pixels, with 
an inner radius of 1.5 x rmax- Any pixels in this area within 
1.5 X Tmax.i of a nearby source, i, or with exposure less than 
10 percent of the median, were rejected. Because of the large 
variations in the PSF, this method works far better than an 
annulus scaled with aperture size and the effect of highly 



varying background, such as around clusters, on the source 
fiux was found to be negligible. 

The source counts are given by 



Counts = Ca — Bkg 

Na Ea 



Bkg = Cb 



Nb Eb 



(1) 
(2) 



where C is the total counts in a region. E is the mean expo- 
sure map value of good pixels, and A'' the number of good 
pixels. Subscripts A and B refer to the source aperture and 
background region respectively. Cs is scaled by the ratio of 
the exposures as in ~ 18 per cent of sources the mean back- 
ground and aperture exposure differ by over 10 percent. 

Throughout the calculations the Gehrels (|Gehrelslll986l ) 
approximation G(C) = 1 + \/0.75 + C is used to approxi- 
mate both the Poissonian 1-sigma upper and lower limit. 
The error on the counts is given by 



0"Co 



{g{Ca)Y + Kkg)= 



(3) 



where the error on the calculated background counts in the 
annulus, (TBkg, is 



2 

0"Bkg 



G{Ce 



Na Ea 
Nb Eb 



+ [G[Cb 



Na Ea 
' Nb Eb 



(4) 



which is the combination of the error on the estimation of 
the background count rate, and the error on applying this 
(low) background value to the aperture. 
The source signal-to-noise ratio (SNR) is then 

SNR = Counts/acounts (5) 

and the significance, SIG (following Ijohnson et al.l (|2003l )) 
is defined as 



so that 



Ca = Bkg + SIG x askg 



SIG = Counts/cTBkg 



(6) 



(7) 



A cut of SIG > 3 was applied to construct a catalogue 
of real sources. A significance of above 3 means that the 
source is not a background fluctuation with above a 3-sigma 
probability. The correlation between SIG and the SNR is 
very good, with a significance cut of 3 corresponding to a 
SNR of around 1.5. This cut is more conservative than the 
WAVDETECT significance parameter, and produces a more 
robust source list. On average it reduces the wavdetect 
source list by around 18 per cent. 

To calculate fluxes (in erg cm~^ s~^, for the 0.5-8 keV 
band) the exposure map value at each pixel (in cm^ s) was 
combined with the counts; 



Flux = — X 
k 



Ed Na Ea st^ Gi 



Nb Eb ^ E, 



(8) 



where the summation over i is over the individual pixels in 
a region, k is the conversion from counts to ergs assuming 
the source has a spectrum with oc v~^''^ between 0.5 and 
8 keV, and e nergy dependent absorption by g alactic hydro- 
gen following Morrison fc McCammon (|l983l ) with column 
density from Dickey fc LockmanI ( 199C ) . 

The fiux missed by choosing a smaller aperture is ~ 1 
per cent for the brightest sources and ~ 4 per cent for the 
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faintest, depending on the source counts. This small cor- 
rection factor was applied to the source fluxes to eliminate 
errors in the full population caused by differences in expo- 
sure times. Luminosities were calculated in the 0.5 - 8 keV 
emission band assuming a Fi, oc v^^''' spectrum. 



5 PREDICTED SOURCE DISTRIBUTIONS 

To interpret the number counts of point sources in each 
image an accurate model of each observation is required 
to determine the number of sources expected if there were 
no AGN in the cluster. This model requires the minimum 
flux detectable at each pixel and the number of blank field 
sources as a function of flux. The changes in sensitivity are 
particularly important in the cluster fields as the extended 
cluster emission may obscure faint central sources. The min- 
imum flux model is described below. Section 15.21 describes 
the calculation of the expected number of sources for each 
observation, and Section [5.31 explains the correction of this 
prediction due to gravitational lensing by the cluster. 

5.1 Modelling the sensitivity of each observation 

A flux limit map wa s computed following the method of 
Ijohnson et al.l (|2003l ). From equation [T] the counts for a 
source centred at pixel i and detected at the minimum signif- 
icance of 3 is C'min,i = 3(TBkg,i, which Combined with equa- 
tion [4] and the conversion to fiux used in equation [8] gives a 
minimum flux detectable with significance > 3 at pixel i of 

„ _ C'min.i 

/ 2 \ 1/2 

= jfc (Rb,.Nb,,) ^) + {GiRB,^NA,,)f] (10) 

where ^min,! is the minimum flux detectable with signifi- 
cance > 3 at pixel i. Subscript A indicates values for the 
predicted source and B the predicted background, and R is 
the rate in counts pixel~^sec~^. The inputs for the predic- 
tion are then the exposure Ei, source size NA,i and back- 
ground count rate Rb.z at each point on the image. 

The exposure is simply the sum of the individual expo- 
sure maps described in Section [JTT] There are regions where 
the gradient in exposure will make source detection difficult, 
and these are masked out later as described below. The er- 
rors in the exposure map should be small and are not easily 
calculable. As they will affect both the blank fields and the 
cluster fields in the same way they can be neglected here. 

For each observation the background rate, including 
the diffuse cluster emission, was calculated by replacing the 
point sources with local background and smoothing the im- 
age with a gaussian kernel of radius 40 pixels. Figure [3] 
shows an example of the background images produced. To 
find the error on the background it is easiest to assume 
that the smoothed background rate is given by the aver- 
age of the counts in a circle, rather than calculating the er- 
rors on the true gaussian convolved image. In other words, 
Rs.i ~ "^ZjizArEa^i/^^^^ wherc the area is a circle of ra- 
dius 40 pixels centred on i. This gives a simple equation for 

the error - aR^ . « -J j . To test the model background. 




Figure 3. Background images for two cluster observations, Abell 
1068 (left) and RX J1720-I-26 (right). The top panel shows the 
the central region of the images, with sources removed and the 
source regions marked by ellipses. The bottom panel shows the 
full smoothed background images (with a square root scale). 



the background rate for each detected source (using aper- 
ture photometry) was compared to that from the smoothed 
images at the same position. The model background accu- 
rately reproduces the calculated background for the detected 
sources (with SIG > 3), with no systematic offset. 

The expected source size distribution was calculated us- 
ing the apertures for the detected sources from 8 represen- 
titive blank fields, and checked against the detected source 
sizes in all fields. Apertures derived from the wavdetect 
output were used instead of the given PSF size as this is 
how the source properties were determined. The radial dis- 
tribution of aperture sizes is shown in Figure [l] which also 
shows the chosen model radial source size distribution. This 
model was determined from the data for significant, low fiux 
sources {S < 0.25 x 10"^'' erg cm~^s~^) which are at the 
detection limit of these observations. The aperture size was 
found to jump at radii of 480, 750 and 1010 pixels, due to the 
behaviour of the PSF combined with the wavelet scales cho- 
sen (this is illustrated by the fact that the brighter sources, 
marked by dots in Figure |31 are far closer to a constant 
slope, as they are affected less by the wavelet scales and 
trace the true change in PSF). The aperture size was mod- 
elled between each jump with a best-fit quadratic, and the 
1-sigma error was determined by the distribution of sources 
around this fit. Above a radius of 1010 pixels the aperture 
size jumps considerably so this area was removed from the 
calculation. 

Again, comparison with the actual source sizes shows 
that this model is accurate to within the errors and has 
no systematic error. There were no significant differences 
between source sizes in the ACIS-I and ACIS-S chips. For 
multiple observations the source size distribution was cal- 
culated for each observation, then combined weighted by 
exposure map. 

A mask was constructed to restrict the area to regions 
where the model is accurate. This removes the effects of 
chip gaps, chip edges and errors in the modelling. Edge ef- 
fects, especially due to the background smoothing, affect 
areas within 60 pixels of the image edge, and 40 pixels of 
chip gaps, so these areas were removed. For merged im- 
ages, i, the 'chip boundary' area was included if ^^(max); < 
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Figure 4. Source size as a function of distance from the opti- 
cal axis (four representitive blank fields are shown to avoid over- 
crowding). The mean size and Ic errors ware calculated for the 
faint sources. This is then the expected source size for sources 
near the flux limit at each pixel. 
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Figure 6. The sky area sensitive to sources of flux > S for one 
field, when no errors are applied to Smin is shown by the dotted 
line and 100 examples of when random errors are applied by the 
solid lines. The variation in the solid lines is negligible compared 
to the error on the number of sources. 




Figure 5. The final flux limit model for Abell 1068 (left) and RX 
J1720-I-26 (right). Light areas are the most sensitive and have 
the lowest limiting flux (1.7 and 1.2 xlO~^^erg cm~^ s~^ re- 
spectively), and dark grey are the least sensitive (7.5 and 1.1 
xlO~^'*erg cm~^ s~^). The PSF size, exposure map and cluster 
background all clearly affect the final limiting flux distribution. 



0.5 X T,j^goodE{max)j , where i5(max) is the on-axis expo- 
sure of an image and the sum is over all images, j, with 
good exposure in the 'chip boundary' area. As described in 
Section 15.11 regions where the model source size is greater 
than 700 pixels were also removed. 

The final flux limit model for the two example flelds is 
shown in Figure [5] where all cuts and masks have been ap- 
plied. To check the flux limit model the fluxes of all detected 
sources were compared to the minimum flux detectable at 
the source position. Almost all (> 97 per cent) of sources are 
brighter than the flux limit at their position. Those that are 
slightly fainter than the corresponding flux limit have large 
errors on their flux, so that <C 1 per cent of sources are over 
la fainter than the calculated flux limit at their position. 

The combined effect of the errors on S'min.i, summed 
over the image, is not straightforward to calculate. Random 
errors were added to the calculation for each pixel, and the 
sky area at each flux re-calculated. As the errors on the 
background level are correlated between pixels, the error in 
each 80 x 80 pixel square were changed by the same (ran- 



domly selected) number of sigma. Changing the size of this 
region did not change the results. The error in calculating 
the flux conversion factor, k, was not included as this will 
affect each fleld in the same way. Figure |S] shows the effect 
of these errors on one field. The sky area without errors has 
quite steep jumps due to the sudden changes in the model 
PSF (due to the wavelet scales) , but applying random errors 
to the flux limit at each pixel smooths this distribution. 

5.2 Log N(>S) — Log S and radial distributions 

The number of sources brighter than a given flux is calcu- 
lated for each field, or for a combination of fields, to produce 
a plot of Log N(>S) against Log S, using 



LogN(> So)=J2f- 

S>So ^ 



(11) 



where N(>So) is the number of sources brighter than So, i 
is the number of sources of flux S, and A the total sky area 
available to detect a source of flux S. 

The errors are dominated by the number of sources de- 
tected, so once the errors have been added to individual pix- 
els (Figure E]), the error on the sky area can be neglected. 
The error on the total number of sources was used for the 
brightest sources, where the sky area is constant, such that; 

o-(is>So) 



0"N(>So) - 



A 



(12) 



for A > 0.99 X A^ax- When the sky area starts to decrease 
(at lower flux), \/i errors are used as they are able to take 
account of the weighting by sky area; at these fluxes the 
number of sources is relatively large (i J>10) and the dif- 
ference between Gehrels and y/i approximations becomes 
minimal. The error is then given by 



CTNOSo) = \ 



S>So 



IS_ 

Ai 



(13) 



It is worth noting here the effect of the Eddington 
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bias l|Eddingtonlll913h . whereby random flux errors can in- 
crease th e meas ured source counts above a chosen flux level. 
iMannersI l|2002h show that the net effect for one fleld is ~ 1 
per cent extra sources above 1.1 x 10~^^ erg cm~^ sec~^. 
Since this is a small effect and will affect the cluster and 
blank fleld samples in the same way, it is not accounted for 
in the analysis. 

For each field, in addition to the Log N(>S) - Log S 
distribution, the radial distribution of sources was found 
and compared to the radial prediction assuming no clus- 
ter sources. This was calculated using the blank field Log 
N(>S) - Log S and the Smin map. Errors on the predicted 
radial distribution were found by applying the Log N(>S) - 
Log S distribution with la errors to the Smin map. The pre- 
dicted and actual radial distributions were calculated from 
the X-ray cluster centre, or from the aim-point if no cluster 
was visible. 

As a check of the pipeline method the radial and Log 
N(>S) - Log S distributions were calculated for the 44 blank 
fields, as described in Appendix A. 3. The pipeline predic- 
tion well reproduces the actual distribution of sources in the 
blank fields. In addition, checks were made for differences 
between the two ACIS detectors on Chandra, as described 
in Appendix A. 4 

5.3 Corrections for gravitational lensing 

The effect of gravitational lensing of X-ray sources by the 
galaxy cluster is small, but is expected to be significant over 
many fields. As discussed by Rcfregicr & Locb ( 1997i ). after 
lensing the fiux of each source is increased by a factor /ig, 
where 6 is the angular distance from the cluster centre, and 
the number density is decreased by the same factor due to 
a decrease in the apparent sky area of the image. Whether 
this results in a net increase or decrease in sources at a 
given flux depends on the slope of the Log N(> 5") - Log 
S distribution. In moderately deep cluster observations the 
slope of the number counts is shallow, resulting in a deflcit 
of sourc es in a lens e d fleld compared to a blank fleld. 

Johnson et all (|2003| ) estimate an expected deficit of 
X-ray sources of ~ 10 per cent in the central 0.5 Mpc of 
MS1054-0321 (z=0.83). This is insignificant for a single field 
but the cumulative effect over many fields may affect the 
sample. In addition, as the effect of lensing on the number 
counts is more significant for bright, moderate redshift clus- 
ters, gravitational lensing could bias the results. 

To calculate exactly the difference between the cluster 
and blank fields that is due to gravitational lensing requires 
detailed knowledge of the dark matter distribution in the 
cluster. As this study is investigating a statistical excess of 
sources in a large number of fields, exact determination of 
the lensing is unnecessary (and unfeasible). Instead, the ra- 
dial loss or gain of sources in each image due to the cluster 
is estimated using a simplified model of gravitational lens- 
ing, with the only inputs being the X-ray luminosity, po- 
sition in the image and redshift of the cluster, an assumed 
background distribution of X-ray sources and the sensitivity 
distribution of the observation. 

For an NFW mass profile (Navarro, Frenk & White, 
1997) is only dependent on a characteristi c radius and 
the cl uster mass, as shown in Append ix A of | Mvers et al.l 
l|2003l ). using formulae and data from iMaoz et aL (^9^ 



iBartelmannI l|l996l 'l and iNavarro et~all (| 19971 '). iMaoz et all 
(|l997l ) also show that the characteristic radius can be ap- 
proximated by a function of the cluster mass. This in turn 
can be estimated by using the redshift de pendent cluster 
mass - luminosity relation in equation 15 of iMaughan et aU 
(2006"). The cluster X-ray luminosities and redshifts (see Sec- 
tions |3]2]T] and |3]2]3]) can therefore be used to calculate the 
distribution of for each field. Although this calculation 
relies on a number of empirical relations, this will not intro- 
duce large errors as discussed below. 

Three models of the X-ray background are used, as de- 
scribed below. They are all calculated for rest-frame 2-8 keV 
(hard band) sources which, where the lensing from clusters 
will be strongest, corresponds to observed 1-4 keV sources. 
The lensing factor calculated in this section is fractional, 
so only the shape and relative normalisation matter. It is 
therefore assumed that the population of hard sources in 
the model shows the same distribution and redshift evolu- 
tion as the sources in the clus ter image. 

Two of the models use the lBarger et all l|2001^ X-ray lu- 
minosity function from 0.1 < z < 1.2, but extend it to z=5. 
The first reduces the density by a factor of at high red- 
shifts, which is a reasonable fit to the sources with confirmed 
redshifts and is therefore a lower limit. The second model 
scales the space density of sources such that the energy den- 
sity per comoving volume remai ns flat a.t z > \. T his is the 
maximum value allowed by thc iBarger et al.l (1200 ll ) data, so 
is an upper limit. The third model adopted here is a luminos- 
ity dependent de nsity evoluti on (LDDE) model, with best fit 
parameters from lUeda et al.l (2003), which is the best fit to 
the hard X-ray luminos ity fu nction from the ChaMP survey 
(see iGreen et all (|2004l ') and [Silverman et all l|2008h for de- 
tails). The three model luminosity functions are calculated 
from z=0 to z=5, in redshift steps of 0.1. The lower end is 
important as a lot of sources will not be lensed, and these will 
reduce any fractional deficit due to lensing. The luminosity 
functions at each redshift are re-normalised to represent the 
sky volume visible in an image of 1 deg^, rather than per 
cubic Mpc. 

The effect of lensing on the model background source 
distribution is calculated as a function of cluster-centric dis- 
tance, cluster redshift and cluster X-ray luminosity for each 
field. The lensed luminosity functions (boosted luminosities 
and lower space densities) of the non-cluster sources were 
found for redshifts - 5 in steps of 0.1 and were converted 
to flux distributions in the observed band and summed over 
all redshifts. The resulting lensed Log N(< S) - Log S dis- 
tributions were compared to the unlensed distribution, and 
the flux limit at each point in the image, to calculate the 
fractional change in sources detected at each pixel. This cor- 
rection was applied up to 300" from the cluster centre, and 
gives a maximum correction per field of ~ 1 source. 

The three models for the X-ray background distribution 
did not give significantly different results (far smaller than 
the errors on the source distribution) so only the LDDE 
model was used, which give s results between t he two ex- 
treme models taken from the lBarger et al.| |200j) data. The 
largest source of error in this model is if there is a system- 
atic miscalculation of the cluster properties, but this is still 
a small source of error overall. For example if the cluster 
mass is assumed to be systematically out by 30 per cent for 
all clusters then this would add the equivalent of 1.5 per- 
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Figure 7. The excess number of sources per field, compared to 
the prediction, for 44 blank fields and 148 uncontaminated cluster 
fields with 0,1 < z < 0.9. The mean for the cluster fields is 1.46 it 
0.32, compared to 0.47 ± 0.61 for the blank fields. 
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Figure 8. The combined Log N(> S) - Log S distribution within 
1 Mpc of the cluster centres for the 44 blank fields and 148 un- 
contaminated cluster fields with 0.1 < z < 0.9. 



cent to the error bar at 1 Mpc. Random errors in the cluster 
properties due to scatter in the cluster scaling relations will 
generally cancel out over a large sample. 

When the correction for gravitational lensing is applied, 
the total number of non-cluster sources predicted in an av- 
erage field decreases. The prediction for the central 3 Mpc 
of the 148 good cluster fields is found to decrease by 0.7 
per cent, which is around 0.27 sources per field on average. 
The calculated number of sources in the cluster, which is the 
number of detected sources minus the prediction, therefore 
increases by the same number. This is shown in Section [6.21 
where the excess sources per cluster field before and after the 
lensing correction are compared. The number of predicted 
sources brighter than 10"^'* erg cm~^sec~^ decreases by 0.6 
per cent, which is 0.06 sources per field. The lensing correc- 
tion is small, typically < 0.5a, but it is not insignificant as 
the correction is predominantly in the central regions. All 
statistics and plots presented in the remainder of this pa- 
per use the lensing correction, but none of the results are 
significantly altered if the lensing correction is ignored. 



6 RESULTS 

6.1 Excess sources in cluster fields 

The number of cluster sources in each field was estimated 
by subtracting the predicted number of sources (Section [5]) 
from the actual number of well-detected sources (Section U]), 
to get the excess sources in each field. The results presented 
in Figure [7] show the excess sources within 1 Mpc of each 
cluster centre, which is the maximum radius observed for 
the lowest redshift clusters. The resulting histogram shows 
that cluster fields have a wide spread of calculated excess 
sources, including negative values, but that the average ex- 
cess is clearly non-zero. For the blank fields, with assigned 
redshifts randomly chosen from the redshift distribution of 
cluster fields, the excess sources within 1 Mpc of each field 
are consistent with zero. Figure [7] shows that the galaxy 
clusters have, on average, around 1.5 sources each within a 



projected distance of 1 Mpc. This value is an average over 
clusters of different redshifts and luminosities, and observa- 
tions of different exposure time; the dependence of number 
of cluster X-ray sources on these variables will be analysed 
in the next paper in this series. 

The Log N(> 5") - Log S distribution was plotted for the 
blank fields and the 148 uncontaminated cluster fields with 
0.1 < z < 0.9. Figure [5] shows that the cluster fields have a 
~ 2a excess at fiuxes of > 10"^'''^ erg cm~'^sec~^. An excess 
of ~ la is seen at fainter fiuxes. These are not strongly de- 
pendent on the lensing correction, which changes the results 
by a maximum of 0.2a". A Kolmogorov-Smirnov (K-S) test 
on the sources brighter than > 10"^"^ '^ erg cm~^sec~^ shows 
that the cluster and blank field populations differ at the 96 
per cent level. 

The most significant excess in the Log N(> S) - Log S 
distribution is found at bright fluxes, but this is partly due to 
the lower number of bright background sources. In fact only 
half of the excess sources within 1 Mpc have flux > lO"^'' 
erg cm~^sec~^ (a mean of 0.76 ±0.18 per field). Sources 
brighter than this fiux are quite likely to be AGN, as this 
corresponds to a (k-corrected) luminosity of > 2.5 x lO**^ 
erg sec~^ in all clusters, and > lO*"^ erg sec~^ in over 80 
per cent of the sample. In clusters with z < 0.2 sources 
with luminosity < 10*^ erg sec~^ can be detected, which are 
far less likely to be AGN, but the majority of sources have 
either fiux ^ 10~^* erg cm^^sec"^ or are in higher redshift 
clusters and so are likely to be AGN. 

6.2 Radial distribution of cluster sources 

The radial distribution of all excess sources in the cluster 
fields is shown in Figure [9j a). The cluster fields clearly have 
an excess of around 1.5 sources per cluster, whereas the 
blank fields have no statistical excess. Although this excess 
is of low significance at 3 Mpc, all of the excess sources lie 
within 1 Mpc from the cluster centres and within this radius 
the significance of the excess is > 3a", with a maximum sig- 
nificance of > 3.5a at 0.85 Mpc. There are no excess sources 
at ^ 1 Mpc despite the fact that two thirds of all detected 
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Figure 9. The number of excess sources per field within a given radius (total detected sources minus the predicted background and 
foreground sources), and la errors, for (a) 148 uncontaminated cluster fields with 0.1 < z < 0.9 corrected for gravitational lousing (red 
line and error bars), and without the lensing correction (green line). All sources are shown, and the cluster fieldss have a significant excess 
of sources within 1 Mpc; (b) the same 148 clusters (thin red line), but restricted to sources with fiux > 10"^'* erg cm^'^sec"^, which 
are detectable in almost all regions of all cluster images. The thick blue 'area corrected' line shows the same distribution corrected for 
the missing area at each radius, due to low sensitivity or gaps in the detector. (The 'area corrected' distribution has higher errors and a 
more variable line at high radius, as a lower fraction of the total area is covered in these regions.); (c) 20 contaminated cluster fields i.e. 
those with another probable cluster in the field of view. As in plot (b), the thin red line shows sources > 10~^* erg cm~^sec~^ and the 
thick blue 'area corrected' line includes a correction for the missing area at high radius. The significant excess of sources at 1 — 3 Mpc, 
when corrected for the missing area, shows that these fields contain sources which are not associated with the main cluster, and therefore 
should not be included in this analysis, (d) The number density of excess sources in each 0.5 Mpc bin, for the same 148 uncontaminated 
cluster fields. The dashed line shows the first bin excluding sources within 25 kpc of the cluster centres. 



sources are beyond this radius. Figure [9ja) also shows the 
radial distribution without correction for gravitational lens- 
ing. As explained in Section [5.31 the number of sources per 
cluster field is lower but still highly significant. 

It is possible that the lack of sources at > 1 Mpc is due 
to the fall off in sensitivity with radius. To check for this the 
radial distribution of sources brighter than 10"^* erg cm~^ 
see" , which can be well detected at all radii, is plotted in 
Figure [9jb). In addition some of the lack of sources could be 
due to the reduced sky area at higher radius, so the distri- 
bution for bright sources was corrected for the proportion of 
missing area at each increase in radius. Both distributions 
show that the cluster sources are still found within ~1 Mpc, 
with a ~ 3a excess in this area. There is no excess above 
this radius, although as the errors are larger in this plot some 



cluster sources could lie beyond 1 Mpc. The correction for 
gravitational lensing for these brighter sources is not plotted 
as it negligible (see Section [5. 3p . 

Figure Oc) shows the same figure for the twenty con- 
taminated fields, which are those which had a second region 
of extended X-ray emission that was not clearly associated 
with the cluster, or an optically detected cluster in the field. 
The distribution of sources is clearly different to that for the 
148 uncontaminated clusters, with excess sources seen up 
to 3 Mpc from the cluster centres. The number of sources 
per cluster in the contaminated fields is similar at 1 Mpc 
(within the errors) but larger at higher radius. This justifies 
the decision not to include these clusters in the analysis, as 
it is very likely that they include sources associated with the 
contaminating clusters on the outskirts of the fields. 
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Figure [9jd) shows the density of sources with flux 
> 10"^'* erg cm~^sec~^ as a function of radius. It is clear 
from this figure that a number of sources lie in the very 
central regions of the cluster, and are likely to be AGN in 
the Brightest Cluster Galaxy (BCG). Wh ereas the fraction 
of radio detected AGN in BCGs is known l|Best et al.ll2007l ') 
the number of X-ray detected AGN is not well defined, and 
detection is complica ted by the extended X-ray emission. 
iRuderman fc Ebelind find that ~10 per cent of their clus- 
ters have a detected X-ray source within 250 kpc of the cD 
galaxy in the cluster centre. 

In this sample 166 clusters have 0.1 < z < 1 and no chip 
boundaries in the cluster centre. Twelve of these 166 clusters 
have X-ray sources within the central 25 kpc, defined from 
the centre of the X-ray emission, whereas only one would be 
expected randomly. One of these clusters, 3C 295, has two 
sources within this radius. Outside this radius there are very 
few additional sources compared to the background predic- 
tion. The k-corrected luminosities of these sources are listed 
in Table (2] To find the proportion of clusters hosting X-ray 
detected AGN it is necessary to find the detection threshold 
at the centre of each cluster. Neglecting the second source 
in 3C 295 the fraction of clusters with AGN with 0.5-8 keV 
k-corrected luminosity > Lx increases from 2A'^\'2 per cent 
at Lx = 10** erg sec~\ to 4.9^1'^ per cent at 10*^ erg sec"^ 
and 7.513 7 P^'" cent at 10*^ erg sec~^. Below this luminos- 
ity most clusters are too bright to detect AGN and so the 
statistics are not significant. 

When the central sources are excluded in Figure [9jd) 
the projected source density is flat or slightly falling until 
~1.25 Mpc, where it falls to zero. This is consistent with 
a random distribution in projected area, which is naively 
not the expected distribution of galaxies in clusters. How- 
ever, the distribution of X-ray sources here is c onsistent with 
the ra dial distribution of cluster galaxies in [Martini et al.l 
l|2007l ). so it may be that X-ray sources simply trace the un- 
derlying population. This will be investigated further in the 
next paper in this series (Gilmour et al. in prep.). 

6.3 Comparison with previous studies 

Most of the clusters that have been studied by other authors 
are also included in this sample, and in the majority of cases 
the results of this analysis agree with the previous results 
within the errors. T he analyses of many small studies, and 
the la rger studies of iMartini et al. (|2007| ) , iBranchesi et al.l 
|2003) and lRuderman fc Ebelind l|2005h . have been repro- 
duced as far as possible and are compared below. 

Comparing the number of excess sources found by 
the pipeline in individual clusters to previous studies, the 
values are in good agreem ent for most clu sters; A2104 
("Martini et al."200^, 3c295 (|Cappi et al."200lV MS1054-03 
( Johnson ct al 2o'o3l'). MS0451-03 (Molna r ct al, 2002 ), six 
of th e clusters fromlCappelluti et all (|2005l ) and MRC 1138- 
262 (|Pentericci et al. 1 120021) . For a few clusters this study 
produces differe nt results in the o ver all number o f sources 
to those found in lCappelluti et al] and lCappi et al.l . as these 
authors present results on a chip-by-chip basis rather than 
as a radial analysis. Visual inspection of these fields indi- 
cates that their conclusions are consistent with those from 
this survey. 

Of the eight clusters studied spectroscopically by 



iMartini et al.l ([20o3), five are included in this study. The 
results agree very well (less than la difference) with the 
IMartini et al.l results, both i n terms of nu mber and radial 
distribution of the sources. IMartini et al.l find 17 sources 
within 1 Mpc in these five clusters, and 12 within 0.5 Mpc. 
This study gives 17.5 within 1 Mpc and 9.5 within 0.5 Mpc. 
At higher radii the errors in this study become too large to 
draw any conclusions for 5 fields. In both studies the AGN 
are predominantly found in the central regions of these 5 
clusters, with twice as many sources at low radius (< 0.5 
Mpc) than at higher radius (0.5 - 1 Mpc). However, from 
Figure IH^b) it is clear that this is not the case for the full 
sample of 148 clusters - rather the number of sources at 
<0.5 Mpc is closer to half the value at high er radius (0.5 
- 1 Mpc). The central concentration of the IMartini et al] 
AGN is therefore no t representativ e of clusters in general, 
perhaps because the IMartini et al.l clusters have particular 
properties. This will be investigated further in the next pa- 
per in this series (Gilmour et al., in prep.) when the 148 
clusters are split into subsamples according to redshift and 
cluster properties. 

IBranchesi et al.l (|2007l ) performed a statistical analysis 
of point sources in 18 clusters, of which 15 are in this sam- 
ple. They find a 1.7(t excess of bright sources within 1 Mpc 
(Fx > 10~^* erg cm~^sec~^ ). The 15 clusters in this pa- 
per have the same 1.7a excess at the bright end of the Log 
N(> S) - Log S distribution, which also fits with th e re- 
sults for the full sample in Figure [8| IBranchesi et al.l find 
7(2) cluster(blank) sources at <0.5 Mpc, and 4(3) sources 
at 0.5-1.0 Mpc and conclude that the majority of sources 
are in the central 0.5 Mpc of the cluster. However, they 
only search to the edge of the intra-cluster emission, which 
gives an average search radius of 0.8 Mpc. The results for 
this study, correcting for missing clusters and different flux 
bands, agree with the Branchcsi ct al] values but the number 
of sources rises steeply beyond 0.8 Mpc, so the conclusion 
that the vast majority of AGN are found within 0.5 Mpc is 
not confirmed if larger radii are investigated. As an aside, 
it is worth noting that five of the fifteen fields investigated 
here are classed as contaminated in this study (morphology 
type Ic, 2c, or 4 in Table [!}. In agreement with Section [6.21 
a significant number of sources continue to be found in these 
fields up to 3 Mpc from the cluster centres. 

iRuderman fc Ebelind (|2005l ) study 51 massive galaxy 
clusters at 0.3 < z < 0.7, and conclude that the point 
sources lie predominantly in the central 0.5 Mpc, with a 
secondary excess at 2 - 3 Mpc. This is significantly different 
from the results shown in Figure HJ^d) . In this study, using 
the 25 clusters with published redshift, the excess in the cen- 
tral 0.5 Mpc is found, but the high si gnificance (8a) of this 
excess found bv lRuderman fc Ebelind and the secondary ex- 
cess at larger radius, a re not. On e possible explanation of 
this is that Rudcrman fc Ebehnd measure their excess from 
the point source density at > 4 Mpc which, as they them- 
selves point out, is lower than that in the control (blank) 
fields. If the value implied from their blank fields is applied 
to the cluster sample, then the secondary excess at 2-3 Mpc 
is no longer significant and the central excess is of lower 
significance, in agreement with this study. It appears that 
their point source density at large radii is artificially low due 
to not subtracting the background sources before scaling to 
physical radius. The density of non-cluster sources, when 
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Figure 10. The 

> 10~^^'^ erg cm~'^sec~^. Solid, hatched and white bars in 
dicate confirmed cluster AGN, possible cluster AGN and back- 
ground/foreground sources respectively. A random distribution, 
based on the available sky area, is marked by the solid line. The 
confirmed and possible cluster AGN are not randomly distributed 
in physical distance from the cluster centres. Ignoring the central 
four AGN, a K-S test shows that the distribution in Mpc is not 
random with 98.5 per cent probability. 



scaled to the cluster redshift, is dependent on the redshift, 
and at higher physical radius only the high redshift fields 
are used to calculate the point source density, leading to a 
lower value. As their point source list and cluster sample are 
not yet published it is not possible to test this further. 



6.4 The brightest sources 

As a further test of the validity of the conclusions the distri- 
bution of the very brightest sources, which are clearly highly 
luminous AGN if they are in the cluster, was investigated in 
detail. Sources brighter than 10~^^'^ erg cm~'^sec~^, which 
are easily detected over the full sample area, were compared 
to the NED to identify possible cluster AGN and eliminate 
contaminating sources. The results, shown in Figure [TOl con- 
firm that the AGN primarily lie within 1.25 Mpc from the 
cluster centre. The 13 AGN which could be cluster members 
are clustered in two groups - one at the cluster centres and 
a second at ~1 Mpc. There is tentative evidence here that 
the very brightest sources lie in the outskirts of the cluster, 
but this will be investigated in the next paper in this se- 
ries. It is clear that these sources are not foreground objects 
but are associated with the clusters, as they are not drawn 
from a random distribution with > 99.93 per cent probabil- 
ity for all AGN, and 98.5 per cent if the central four AGN 
are excluded. 



7 CONCLUSIONS 

The X-ray point source population in moderately deep (> 10 
ksec) Chandra observations of 148 cluster fields and 44 blank 
fields were calculated and compared in order to estimate the 
number of X-ray sources in the galaxy clusters. The number 
of sources is found to be low, with 1.5 sources in a typical 



cluster. This result is significant to > 3cr, but the actual 
number of sources per cluster depends on the exposure time 
and the redshift of the cluster. Over half of these sources 
have fiuxes corresponding to luminosities > 10*'^ erg sec~^ 
and are likely to be AGN. The population is not dominated 
by AGN in the central galaxy (BCG), as only twelve clus- 
ters have a central source, rather the sources are AGN or 
star-burst galaxies in normal cluster members. They are all 
found in the central 1 Mpc of the cluster, and are randomly 
distributed in projected area within this radius. 

Many of the cluster s cove re d by similar stud- 



ies (e.g. iBranchesi et all l2007l : iMartini etahl l2007l : 



iRuderman fc Ebehnd |2005| ) are included in this sam- 



ple, and when the same clusters are compared then the 
results are generally in good agreement. However the 
conclusions drawn from these papers, which use smaller 
numbers of clusters, are often not borne out by this larger 
study. For example a higher number of sources is found in 
the central 0.5 Mpc than the annulus at 0.5 - 1 Mpc in 
all three previous studies of more than six clusters, and 
whilst these results can be reproduced by this study for the 
smaller samples, the cluster population in general has more 
sources in the outskirts (0.5 - 1 Mpc) of the cluster than the 
central 0.5 Mpc. In Appendix A. 2 it is demonstrated that 
samples of less than five clusters suffer strongly from cosmic 
variance in the number of background sources, but larger 
samples are also affected to some extent. The question of 
whether the discrepancies between this study and some of 
the previous papers in this field is due to cluster properties 
or the larger sample size in this study will be investigated 
in the next paper in this series. 

This paper describes a sample of point sources in clus- 
ter fields which can be used to investigate the number and 
properties of X-ray sources in galaxy clusters as a function 
of cluster properties and redshift, and hence increase our 
understanding of the links between environment and AGN. 
This, the first paper in the series, serves as an introduction 
to the sample and comparison to previous studies. Subse- 
quent papers (Gilmour et al. in prep) will use this study to 
investigate in more detail the environments of cluster AGN. 
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Table 2 



NED Name Obs Ids RA DEC Array Exp L^l L^2 rx^ z Ref Excess Centre 



Morphology 1 clusters: 
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00 
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.6 
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.2 


2 
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23 
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8 


.31 
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1,2 


-2.37+j- 
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00 
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.7 


-12:22:41.9 


1 


24, 


.7 


1 
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13 


.40 


.5 


79 
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3 
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00 
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35 


.8 
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S 


38, 


.9 
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2 
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1 
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90 
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:27: 
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.4 
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I 


9, 


.9 


0, 
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1 
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0. 
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5 
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;34, 


.0 
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0. 


.94 
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I 
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.0 
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0. 
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7 
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.4 
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I 


19 


.9 
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:59; 
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.1 
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I 
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.5 
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36 


.1 
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1. 
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9 


6.47lj; 


CL J023026.6+183622 


5754 


02; 


:30; 


28, 


.6 


+ 18:36:14.7 


I 


67. 


.6 


0, 


.14 


3, 


.88 


1. 


.36 
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.8 
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.6 
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3 
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1 
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2 
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.7 
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.0 


1 
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8, 
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1 
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1 
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07 
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.4 
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s 
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.8 
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14 
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4, 
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1 
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19 
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:09: 
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.9 


+28:12:06.9 


I 
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2 
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1 
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0, 
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2 
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.0 
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S 
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.5 


2 
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1 
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1 
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16 
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.0 
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s 


8 


.2 


1 


.63 
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.33 


1 


.02 
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20 
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08 


:30; 


;58 


.7 
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1 


38, 


.6 


9 


.12 


10 
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4, 


.50 


1 


.13 


0.182 


7,21 
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4C +55.16 
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08 


:34; 


,5,5 


.0 


+ 55:34:21.2 


6,7 


92, 


.0 


3 


.59 


6 


.35 


3, 


.29 


1 


.44 


0.242 


22 
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2MASX .108425596+2927272 
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08 


:42- 


55 


9 


+29:27:25.5 


S 


29, 


.4 


3 


.14 


3 


.39 


1 


.84 


1 


.25 
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2 


-3.56l^; 
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08 


:42 


57 


.5 


+36:21:56.1 


1 


23, 


.6 


7 


.30 


16 


.61 


6, 


.50 


1 


.10 


0.282 


7 
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RX J0850.1+3604 


1659 


08 


:50; 


;06 


.7 


+ 36:04:17.1 


1 


22 


.1 


2 


.94 


12 


.99 


5 


.19 


0, 


.93 


0.378 


5,16 


2.681*; 


ZwCl 0848.5+3341 


4205 


08; 


:51: 


39, 


.0 


+33:31:08.0 


S 


11. 


.4 


1 


.03 


4, 


.81 


2, 


.04 


1 


.15 


0.380 


23 


-2.22l^- 


MACS J0913.7+4056 


509 


09; 


:13; 


;45, 


.4 


+40:56:27.6 


s 


7, 


.8 


1 


.46 


10, 


.00 


4, 


,18 


0. 


.87 


0.442 


24 


0.213- 


ABELL 0773 


533,3588,5006 


09; 


:17; 


:52, 


.9 


+51:43:39.2 


I 


30. 


.4 


6, 


.57 


8, 


.54 


3, 


.87 


1 


.03 


0.217 


7,5 


10.281®; 


RX J0926.6+1242 


5838 


09; 


;26; 


36, 


.6 


+ 12:43:03.4 


I 


31. 


.3 


0, 


.30 


2, 


.77 


1 


.41 


1 


.02 


0.489 


8 


3.921*; 


MACS J0947. 2+7623 


2202 


09; 


;47; 


13, 


.0 


+76:23:14.2 


I 


11. 


.7 


6, 


.30 


24, 


.61 


11. 


.48 


0, 


.91 


0.36 


25 


1.1513; 


ZwCl 0947.2+1723 


3274 


09; 


;49; 


;51, 


.8 


+ 17:07:08.1 


I 


14. 


.3 


3, 


.19 


14, 


.28 


5, 


.42 


0, 


.73 


0.383 


5 


3.9813; 


ZwCl 0949.6+6207 


3195 


09; 


;52; 


49, 


.3 


+51:53:04.9 


s 


26. 


.6 


4, 


.35 


5, 


.87 


3, 


.26 


1 


.01 


0.214 


16 


6.88l|; 


RX J0956.0+4107 


5759 


09; 


:56; 


03, 


.2 


+41:07:13.0 


I 


40. 


.0 


0, 


.28 


3 


.79 


1 


.81 


0, 


.98 


0.587 


8 


-A.2+\ 


CL J095819. 3+470217 


5779 


09; 


:58; 


19, 


.2 


+47:02:03.5 


I 


26. 


.1 


0, 


.25 


1 


.40 


0, 


.82 


1 


.06 


0.39 


26 


0.8913; 


ABELL 0907 


535,3185,3205 


09; 


:58: 


:21 


.9 


-11:03:50.9 


I 


104, 


.6 


8 


.58 


5 


.31 


2 


.44 


1 


.42 


0.153 


27 


7.60l^; 


MS 1008.1-1224 


926 


10; 


:10; 


:32, 


.3 


-12:39:34.5 


I 


43, 


.6 


2, 


.21 


6, 


.00 


2, 


.40 


1 


.23 


0.301 


2 


3-931*; 


ZwCl 1021.0+0426 


909 


10; 


:23: 


39, 


.7 


+04:11:09.2 


I 


46. 


.8 


10. 


.14 


25 


.66 


11. 


.40 


1 


.00 


0.291 


28 


4.78l=- 


ABELL 1068 


1652 


10 


:40; 


44 


.6 


+39:57:10.2 


s 


26, 


.7 


10, 


.63 




.41 


3, 


.22 


1 


.50 


0.138 


7 


8.60I?/ 


RX J1008. 8+0906 


3252,5009 


11 


:08; 




.3 


+ 09:05:58.5 


1 


34, 


.1 


1 


.18 


8 


.73 


3, 


.77 


0, 


.86 


0.463 


23 


0.46+3- 


WARP .11113.0-2615 


915 


11 


:13; 


05 


.1 


-26:15:38.8 


1 


104, 


.2 





.10 


2 


.46 


1 


.01 


1 


.09 


0.725 


9 


i.9i*: 


ABELL 1204 


2205 


11 


:13 


20 


.4 


+ 17:35:38.7 


I 


23 




6 


.78 


.5 


.52 


3, 


.41 


1 


.24 


0.171 


7,5 


i.89l|; 


RX .11115.8+0129 


3275 


11 


:15; 


51 


.9 


+ 01:29:55.6 


1 


14, 


6 


4 


.84 


18 


.57 


7, 


.72 


0, 


.93 


0.38 


29 (Nl) 


0.9513; 


RX ,11120.1+4318 


5771 


11 


:20; 


06 


.9 


+43:18:06.5 


1 


19 


8 





.41 


6 


.11 


3, 


.00 


0, 


.87 


0.6 


30 


1.31+3- 


RX J113(). 9+2326 


1660 


11 


:20; 


57 


.4 


+23:26:33.1 


1 


70 


6 





.25 


3 


.30 


1 


.79 


0, 


.95 


0.562 


8 




MS 1137.5+6624 


536 


11 


:40: 


22, 


.3 


+ 66:08:16.1 


I 


117. 


.2 


0, 


.27 


7, 


.21 


2, 


.94 


0, 


.87 


0.782 


31 


-3.48+2- 


ABELL 1361 


2200,3369 


11; 


:43; 


:39, 


.7 


+46:21:20.0 


s 


16. 


.9 


4, 


.60 


1 


.64 


0, 


.97 


1 


.07 


0.117 


7,5 


0.8l'l- 


ABELL 1413 


537,1661,5003 


11; 


:55; 


:18, 


.0 


+23:24:16.2 


I 


94 


.1 


16. 


.17 


8, 


.33 


3 


.65 


1. 


.28 


0.142 


7,5 


o.55l«; 


ABELL 1446 


4975 


12; 


:02; 


:04, 


.7 


+58:02:12.0 


s 


58, 


.2 


2, 


.87 


0, 


.77 


0. 


.44 


1. 


.11 


0.103 


5 


o.36l^; 


CLG J1205+4429 


4162 


12; 


:05; 


;51, 


.4 


+44:29:10.8 


s 


29. 


.7 


0, 


.05 


0, 


.82 


0, 


.50 


1 


.08 


0.592 


32 


-3.88+j- 


RXC J1206. 2-0848 


3277 


12; 


:06; 


12, 


.4 


-08:48:03.9 


I 


23. 


.4 


4, 


.70 


28, 


.82 


10. 


.65 


1 


.07 


0.441 


33,29 


1.4ll3i 



59 
.50 

ie 41.70 

94 
.66 
60 
.40 

Continued on next page. . 
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Table 2 - Continued 



NED Name Obs Ids RA DEC Array Exp Fj, L^l L^2 rx^ z Ref Excess Centre 



RBS 1080 


5833 


12 


13 


23 


1 


-26:18:07.9 


6,7 


9 


9 





61 


1 


41 





60 





96 


0.278 


34 




40 


43.16 


RX J1213. 5+0253 


4934 


12 


13 


35 





+02:53:47.9 


I 


18 


7 





25 


1 


38 





63 





77 


0.409 


8 


0.60+^ 


60 




R,X .11216.3+2633 


4931 


12 


16 


19 


9 


+ 26:33:12.4 


I 


17 


5 





20 


1 


41 





80 





86 


0.42S 


6 




18 




RX ,71221.4+4918 


1662 


12 


21 


26 


3 


+49:18:27.2 


I 


78 


7 





38 


7 


47 


3 


26 





95 


0.7 


6 


0.41 + ^ 


78 




CL .7122201.9 + 270919 


5766 


12 


22 


01 


9 


+ 27:09:32.4 


I 


49 








18 


1 


54 





85 


1 


OS 


0.472 


10 


0.49j:^ 


96 




BMW-HRI ,7122657.3 + 333253 


5014 


12 


26 


58 





+ 33:32:47.1 


I 


32 


6 





86 


26 


81 


10 


23 


1 


04 


0.S9 


35 


1.23+1 


39 




ABELL 3541 


1648 


13 


03 


42 


4 


-24:14:45.3 




9 


7 


7 


78 


3 


35 


1 


80 


1 


28 


0.128 


7 


-4.45+^ 


SO 




MACS ,71311.0-0311 


3258,6110 


13 


11 


01 


7 


-03:10:38.5 


I 


77 


8 


1 


11 


9 


92 


4 


79 





95 


0.49 


13 


4.6ll* 


71 




ABELL 1689 


7289 


13 


11 


29 


5 


-01:20:29.8 


I 


74 


9 


18 


52 


16 


40 


7 


10 


1 


42 


0.183 


7 


5.92+^ 


57 




RX .71320.0+7003 


3278 


13 


20 


07 


9 


+70:04:36.8 


I 


20 


5 


1 


74 


5 


74 


2 


58 





73 


0.328 


7,5 


5.42l^ 


57 




ZwCl 1332.8+5043 


5772 


13 


34 


20 


1 


+ 50:31:01.2 


I 


17 


6 





36 


5 


54 


2 


67 


1 


00 


0.62 


30 


3.7ll^ 


78 




RX J1340.5+4017 


3223 


13 


40 


32 


9 


+40:17:38.7 




46 


4 





16 





13 





13 


1 


06 


0.171 


36 


6.78+1 


18 




LCDCS 0829 


3592 


13 


47 


30 


8 


-11:46:10.1 


I 


67 


4 


10 


78 


68 


09 


23 


48 


1 


16 


0.461 


37,38 


-1.92+^ 


40 




RDCS J1350+6007 


2229 


13 


50 


48 


3 


+60:07:06.0 


I 


68 


1 





12 


3 


99 


1 


98 


1 


06 


0.804 


39 


2.70j:2 


12 




ZwCl 1358.1+6245 


516 


13 


59 


50 


6 


+62:31:02.9 


S 


62 


1 


2 


42 


8 


31 


3 


80 


1 


11 


0.328 


2 


1.52+^ 


34 




ABELL 1835 


495,496 


14 


01 


02 





+02:52:41.6 


S 


30 


2 


18 


40 


34 


36 


14 


89 


1 


68 


0.263 


7 


-2.33+1 


if 




3C 295 


578 


14 


11 


20 


4 


+ 52:12:10.0 


s 


17 


9 





87 


6 


92 


3 


43 


1 


21 


0.46 


40 


r.57tt 


84 


43.77(N2) 


NSCS J141623+444558 


541 


14 


16 


27 


9 


+44:46:44.5 


I 


30 


7 





58 


3 


38 


1 


80 


1 


02 


0.4 


6 


-2.53^2 


40 




MACS J1423. 8+2404 


4195 


14 


23 


47 


9 


+24:04:42.6 


6,7 


116 


2 


1 


92 


22 


14 


9 


91 


1 


20 


0.646 


1 


3.73j:| 


10 




ABELL 1914 


542,3593 


14 


26 


02 





+37:49:32.8 


I 


26 


7 


19 


03 


14 


54 


6 


18 


1 


00 


0.171 


7 


6.23+1 


78 




RBS 1460 


5793 


15 


04 


07 


5 


-02:48:16.1 


I 


38 


9 


25 


15 


32 


46 


13 


57 


1 


22 


0.215 


41 


-lists 


72 




ABELL 2034 


2204 


15 


10 


11 


8 


+33:30:54.3 


1 


53 


8 


10 


73 


3 


41 


1 


59 





97 


0.113 


7 


3.82^7 


49 




RX J1532. 5+3021 


1665 


15 


32 


53 


8 


+30:20:58.6 


7 


9 


9 


,5 


92 


22 


92 


11 


18 


1 


OS 


0.345 


16 (N3) 


1.18+2 


40 




ABELL 2111 


544 


15 


39 


41 


3 


+34:25:06.7 


7 


10 


3 


3 


57 


5 


18 


2 


26 


1 


02 


0.229 


7 


-1.70+1 


40 




ABELL 2104 


895 


15 


40 


07 


9 


-03:18:17.5 


S 


49 





8 


94 


5 


42 


1 


80 


1 


72 


0.153 


7,42 


12.18+g 


41 




WARP ,71552.2+2013 


3214 


15 


52 


12 


8 


+20:13:39.9 


S 


14 


9 





23 





11 





08 


1 


06 


0.136 


6 


-3.30tJ 


It 




MACS ,71621.4+3810 


3594,6109,6172 


16 


21 


24 


7 


+ 38:10:08.5 


7 


77 


2 


1 


50 


11 


50 


5 


45 


1 


03 


0.465 


43 


0.12+^ 


96 




ABELL 2204 


499 


16 


32 


46 


9 


+05:34:31.5 


S 


10 





32 


58 


19 


81 


8 


68 


1 


24 


0.152 


7 


1.89+1 


58 




ABELL 2218 


553,1454,1666 


16 


35 


51 


5 


+ 66:12:36.8 


7 


60 


3 


5 


76 


4 


77 


2 


11 





96 


0.176 


7 


7.1+1 


92 




ABELL 2219 


896 


16 


40 


19 


9 


+46:42:35.3 


S 


42 


1 


16 


04 


22 


61 


9 


30 


1 


57 


0.226 


7 


-7.r>itt 


30 




Hercules A 


6257 


16 


51 


08 


2 


+04:59:33.0 


6,7 


49 


3 


4 


34 


2 


77 


1 


40 


1 


16 


0.154 


44 


2.96+'3 


7.5 




ABELL 2256 


3245 


17 


20 


08 


4 


+ 27:40:11.0 


I 


10 





6 


37 


4 


58 


2 


14 





99 


0.164 


7 


2.35 + 1 


44 




SDSS-C4 3072 


1453,3224,4361 


17 


20 


10 


1 


+26:37:30.8 


7 


64 


4 


14 


37 


10 


30 


4 


69 


1 


22 


0.164 


16 


1.2+1 


20 




MACS J1720. 2+3536 


6107 


17 


20 


16 


7 


+35:36:23.9 


7 


33 


6 


2 


59 


13 


05 


6 


69 


1 


07 


0.391 


43 


0.93+^ 


96 




ABELL 2294 


3246 


17 


24 


12 





+85:53:10.8 


7 


9 


6 


7 


97 


6 


48 


2 


23 


1 


06 


0.178 


7 


9.13+1 


22 




MS 2053.7-0449 


551,1667 


20 


56 


21 


1 


-04:37:46.3 


7 


88 


4 





21 


2 


72 


1 


13 





96 


0.683 


2 


1.83+3 


43 




MACS J2129. 4-0741 


3199,3595 


21 


29 


26 





-07:41:28.2 


7 


36 


7 


1 


61 


18 


25 


6 


67 


1 


00 


0.67 


1 


6.77iJ 


57 




RBS 1748 


552 


21 


29 


40 





+00:05:19.7 


7 


9 


9 


6 


72 


9 


54 


4 


24 





98 


0.224 


16 


0.70+^ 


79 

60 




MS 2137.3-2353 


928 


21 


40 


15 


2 


-23:39:40.1 


S 


39 


4 


3 


69 


11 


63 


6 


64 


1 


34 


0.313 


2 


-4.44t^ 


44 
30 




ABELL 2390 


4193 


21 


53 


39 





+ 17:41:15.3 


6,7 


93 


7 


16 


19 


23 


69 


9 


01 


1 


66 


0.23 


7,45 


12.22+5 




92 


42.90 


ABELL 2409 


3247 


22 


00 


52 


9 


+20:58:22.3 


7 


10 


2 


8 


75 


5 


00 


2 


15 





97 


0.148 


7 


-3.44+^ 


41 
20 




IRXS J221144. 6-034947 


3284 


22 


11 


45 


9 


-03:49:46.6 


7 


17 


7 


7 


06 


14 


29 


5 


13 


1 


04 


0.27 


41 


7.5tt 


84 
71 




MACS J2214. 9-1359 


3259 


22 


14 


57 


3 


-14:00:12.3 


7 


19 


2 


1 


92 


15 


07 




98 


1 


09 


0.483 


46 


2.81^2 


78 
60 




MACS .72229.8-2756 


3286 


22 


29 


45 


2 


-27-55-36.4 


7 


16 





3 


22 


11 


12 


6 


36 





90 


0.322 


47 


2.13+^ 


96 
78 




IRXS ,7224322.6-093549 


3260 


22 


43 


20 


9 


-09:35:42.8 


7 


20 


3 


2 


72 


17 


42 


7 


02 


1 


06 


0.439 


48 


1.52^1 


60 
40 




IRXS ,7224505.2+263758 


3287 


22 


45 


04 


7 


+ 26:38:03.5 


7 


16 


1 


3 


36 


9 


52 


4 


02 


1 


05 


0.304 


,5 


0.83^2 


40 
18 




RX ,72247.2+0337 


911 


22 


47 


28 





+ 03:37:00.8 


7 


48 


7 





12 





11 





08 


1 


02 


0.18 


6 


14.40+g 


48 
41 




ABELL S1063 


4966 


22 


48 


44 


8 


-44:31:46.4 


7 


26 


6 


11 


12 


40 


28 


16 


32 


1 


09 


0.348 


29 


7.,30+* 


71 
57 




RBS 1906 


5769 


22 


51 


47 


,5 


-32:06:12.5 


6,7 


10 


3 





32 





53 





21 





74 


0.246 


41 


7.40+= 


9 

97 


44.00 


ABELL S107 


1562 


22 


58 


48 


2 


-34:48:07.3 


S 


72 


1 


4 


74 


13 


72 


5 


69 


2 


18 


0.31 


49,50 


2.70t| 


68 
57 




ABELL 2537 


4962 


23 


08 


22 


1 


-02:11:27.9 


6,7 


36 





2 


85 


7 


27 


2 


88 


1 


19 


0.296 


61 


-3.29+* 


44 

29 


42.04 


ABELL 2631 


3248 


23 


37 


38 


7 


+ 00:16:08.5 


7 


9 


1 


4 


12 


9 


19 


3 


78 


1 


05 


0.278 


7.29 


4,65+* 


12 




ADELL 2(.(.7 


2214 




."li 






-2(j:U,"'kU:',,,") 




!) 


(> 


li 


!)IJ 


17 




^ 


71 


1 


14 


IJ.226 


52.2!) 




(ii 




Morphology 2 clusters: 


RX J0404.6+1109 


3269 


04 


04 


32 


9 


+ 11:08:08.0 




21 


7 





89 


3 


42 


1 


14 





96 


0.366 


5 


5.69+* 


43 
28 




RX J0853.2+5759 


5765 


08 


53 


16 


8 


+57:59:44.4 




24 


6 





16 


1 


40 





72 





80 


0.476 


6 


2.73+1 


96 
78 




RX J1006.9+3200 


5819 


10 


06 


54 


5 


+32:01:32.6 




10 


8 


1 


45 


7 


58 


3 


61 


1 


08 


0.398 


64 


1.68j:^ 


39 
18 




ZwCl 1006.1+1201 


925 


10 


08 


47 


5 


+ 11:47:36.1 




29 


3 


2 


63 


3 


61 


1 


61 


1 


26 


0.221 


2 


7.61+1 


56 
45 




ABELL 1201 


4216 


11 


12 


54 


6 


+ 13:26:02.4 




34 


2 


3 


21 


2 


51 


1 


29 


1 


22 


0.169 


7 




24 
13 




ZwCl 1112.2+5318 


5008 


11 


15 


15 


8 


+53:19:54.3 




18 





1 


56 


11 


65 


6 


60 





87 


0.466 


23 


1.7615 


60 
40 




ABELL 1240 


4961 


11 


23 


37 


3 


+43:06:08.3 




61 


1 





93 


1 


01 





64 





97 


0.196 


66 


-3.51+1 


80 
70 




NSCS J125606+255746 


3212 


12 


56 


02 


5 


+25:56:38.0 




26 


9 





13 





22 





14 





91 


0.232 


6,66 


-6.3li*3 


59 
45 




ABELL 1763 


3591 


13 


35 


18 


3 


+41:00:00.7 




19 


6 


5 


85 


8 


08 


3 


78 


1 


26 


0.223 


67 


-1.78+1 


96 
79 




RX J1354.2-0221 


5835 


13 


54 


17 


1 


-02:21:52.6 




37 


6 





14 


1 


85 


1 


01 





87 


0.546 


8 


-1.13tl 


40 

18 




MS 1621.5+2640 


546 


16 


23 


35 


2 


+26:34:21.3 




29 


9 


1 


10 


6 


54 


2 


60 


1 


02 


0.426 


2 


0.94+5 


96 
78 




RX J1716. 4+6708 


548 


17 


16 


48 


9 


+67:08:25.4 




51 


6 





32 


8 


82 


3 


18 





76 


0.813 


58 


^ 4.491* 


28 
12 
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Tabic 2 - Continued 



NED Name Obs Ids RA DEC Array Exp L^l L^2 rx^ z Rof Excess Centre 



ABELL 2261 


550,5007 


17:22:27.2 


+32:07:57.0 


I 


33.3 


9.02 


12.53 


5 


16 


1.08 


0.224 


7 


MACS J1824. 2+4309 


3255 


18:24:18.5 


+43:09:54.2 


I 


14.9 


0.03 


0.41 


0.36 


0.07 


0.487 


43 


MACS J2228. 5+2036 


3285 


22:28:33.2 


+20:37:12.9 


I 


19 


8 


2.99 


16.25 


6.30 





81 


0.412 


5 


ABELL 2550 


2225 


23:11:35.7 


-21:44:46.8 


s 


58 


6 


1.03 


0.42 


0.34 


1.23 


0.123 


59 


Morphology 3 clusters: 


ABELL 2744 


2212 


00:14:19.1 


-30:23:23.1 


6,7 


24 


7 


8.18 


23 


31 


9 


55 


1 


25 


0.308 


7 


ABELL 0115 


3233 


00:55:50.5 


+26:24:35.5 




49 


6 


2.98 


3 


28 


1 


65 


1 


02 


0.197 


7 


EDCC 586 


5778 


01:41:32.6 


-30:34:42.6 




29 


3 


0.17 





14 





13 


1 


23 


0.168 


60 


CL J0152. 7-1357 


913 


01:52:44.4 


-13:57:17.4 




36 


3 


0.19 


5 


12 


1 


94 





99 


0.831 


61 


ABELL 0521 


901 


04:54:06.6 


-10:13:09.4 




38 


5 


3.50 


6 


48 


2 


72 


1 


09 


0.253 


62 


ABELL 0520 


528,4215 


04:54:10.1 


+02:54:40.7 




75 


5 


6.38 


7 


14 


2 


81 


1 


20 


0.203 


7,5 


IRXS J065830. 3-555702 


554,3184 


06:58:30.3 


-55:56:34.6 




112 


4 


13.22 


32 


58 


11 


24 


1 


14 


0.296 


63 


CXOU J091554+293316 


4209 


09:15:52.5 


+29:33:23.6 




19 


1 


0.08 





83 





42 





85 


0.500 


64 


RXC J1234. 2+0947 


539 


12:34:21.7 


+09:46:56.9 




9 


1 


2.43 


3 


66 


1 


83 





79 


0.229 


5 


ABELL 1682 


3244 


13:06:50.5 


+46:33:25.3 




8 


9 


4.51 


6 


83 


2 


86 


1 


12 


0.234 


7 


ABELL 1758 


2213 


13:32:42.8 


+50:32:55.0 


6,7 


55 


8 


7.49 


17 


28 


7 


26 


1 


43 


0.28 


7,5 


ABELL 2069 


4965 


15:24:08.7 


+29:53:00.3 




52 


2 


3.19 


1 


09 





56 


1 


08 


0.116 


7,65 


FIRST .1234229.5+001845 


5786 


23:43:41.9 


+ 00:18:07.5 




29 


7 


1.49 


3 


22 


1 


42 


1 


08 


0.27 


13 



-6.87l3;61 
0.3913.18 

O c:o + 3-78 
'=-O^_2.60 

-4 40 + '' '''* 
^•*"-5.70 



5 96+''-^'' 

1 Q 1 ~l" 3 ■ 3^ 

l.ol_2 
., „, + 6.19 

' ■ '"-4,56 
, + 3.39 

0.77tl- 

"^-4.11 
-4.821^-3 



Morphology Ic clusters: 



ABELL 0267 


1448,3580 


01:52:42.2 


+ 01:00:40.5 


27 


4 


4 


16 


6 


27 


2 


81 


1 


13 


0.23 


7,5 


3.14+^ 


TFT — 

85 




MACS J0159. 8-0849 


3265,6106 


01:59:49.3 


-08:49:59.9 


52 


8 


4 


18 


21 


81 


9 


66 


1 


06 


0.4 


48 


2.161* 


12 
96 




RDCS 0337.4-3457 


6264 


03:37:25.0 


-34:57:18.6 


12 


2 





02 


1 


30 


1 


06 





26 


0.84 


66 


0.10+^ 


65 
32 




CL J034051. 6-282310 


5780 


03:40:52.9 


-28:23:08.3 


24 


6 





37 


1 


53 





94 


1 


02 


0.346 


10 


2.2tl 


28 
12 




RX J0451. 9+0006 


5815 


04:51:54.4 


+ 00:06:19.2 


10 


2 


1 


36 


8 


27 


3 


05 





81 


0.43 


67 


3.34^^ 


59 
39 




MACS J0454. 1-0300 


902 


04:54:11.2 


-03:00:51.3 £ 


43 


5 


2 


37 


25 


07 


9 


06 


1 


08 


0.55 


1,31 


-4.82^^ 


94 
66 




MACS J0717+3745 


1655,4200 


07:17:31.3 


+37:45:29.5 


78 


7 


3 


69 


35 


91 


11 


53 


1 


17 


0.548 


1 


3. 80+* 


28 
12 




ABELL 0586 


530 


07:32:20.3 


+31:37:56.2 


10 





7 


14 


5 


51 


2 


27 


1 


02 


0.171 


7 


0.96t* 


13 
97 




ABELL 0611 


3194 


08:00:56.7 


+36:03:23.2 £ 


35 


9 


3 


82 


9 


62 


4 


18 


1 


22 


0.288 


7 


2.99^4 


22 
10 


43.29 


MS 0906.5+1110 


924 


09:09:12.8 


+ 10:58:33.0 


29 


6 


4 


69 


3 


85 


1 


73 





94 


0.175 


5,2,7 


1.40t= 


22 
11 




ABELL 0963 


903 


10:17:03.6 


+39:02:53.4 £ 


36 


2 


7 


98 


9 


54 


4 


53 


1 


16 


0.206 


7,5 


7.75+^ 


38 
29 




MS 1054-03 


512 


10:56:59.0 


-03:37:35.0 £ 


85 


7 





63 


18 


13 


6 


77 


1 


04 


0.823 


31,68 


0.42t* 


13 
96 




MACS J1149. 5+2223 


1656,3589 


11:49:35.0 


+22:24:06.7 


38 


5 


2 


20 


1 


82 





83 


1 


06 


0.176 


5,1 




19 
10 




CL J131219. 4+390058 


5781 


13:12:19.5 


+39:00:53.4 


25 


1 





23 


1 


46 





93 


1 


30 


0.404 


10 


-3.52+^ 


94 
66 




RDCS J1317+2911 


2228 


13:17:21.5 


+ 29:11:16.2 


111 








01 





58 





31 





72 


0.805 


39 


4.2+* 


71 

58 




MS 1455.0+2232 


543,4192 


14:57:15.0 


+22:20:34.7 


101 


4 


5 


97 


11 


96 


6 


04 


1 


28 


0.258 


28 


-5.33t* 


73 

60 




ABELL 2163 


545,1653 


16:15:45.8 


-06:09:01.8 


80 


3 


24 


45 


25 


75 


7 


58 


1 


65 


0.203 


69,70 


-6.78+* 


30 
15 




IRXS J201127.9-572507 


4995 


20:11:27.1 


-57:25:10.1 


23 


9 


1 


55 


3 


75 


1 


91 





77 


0.279 


41 


1.59t* 


28 
13 




Morphology 2c clusters: 


ABELL 1300 


3276 


11:31:54.8 


-19:55:49.5 


13.8 


4.47 


12.36 


4.74 


0.90 


0.307 


7,71 


„ , +3.9B 
•^■.^.^-2.78 




ABELL 1942 


3290 


14:38:21.9 


+03:40:09.7 


56.7 


1 


14 


1.56 


0.73 


0.96 


0.22 


72 


7Atl 


18 
9 


41.93 



High redshift (proto)clusters: 



XLSSC 029 


7185 


02 


24 


04 


1 


-04:13:30.1 


6,7 


32 


9 






— — 1.05 


73 


2 


20+2 


78 

60 


43.24 


RCS J0439-2904 


3S77 


04 


39 


38 


1 


-29:04:55.8 


S 


86 


6 






— — 0.951 


53 





19+^ 


60 
40 




CL J0442+0202 


3242 


04 


42 


23 


8 


+02:02:19.6 


I 


43 


4 






— — 1.11 


74 


6 


15+* 


27 

12 


44.75 


3C 184 


3226 


07 


39 


28 


4 


+70:23:39.8 


S 


18 


8 






— — 0.996 


75 





67+1 


18 
93 




3C 210 


5821 


08 


58 


10 





+27:50:53.5 


6,7 


20 


5 






— — 1.169 


76 





50+1 


93 
66 


44.06 


RDCS J0910+5422 


2227,2452 


09 


10 


44 


6 


+54:22:03.8 


I 


168 


8 






— — 1.11 


77 


9 


2i+_\ 


22 
10 




CL J100207. 7+685848 


5773 


10 


02 


09 


2 


+68:58:38.0 


I 


19 


8 






— — 0.928 


10 


4 


41+_l 


78 
59 




PKS 1138-26 


898 


11 


40 


48 


3 


-26:29:09.9 


S 


32 


5 






— — 2.16 


78 


4 


83+* 


28 
12 


45.58 


RDCS J1252-2927 


4198,4403 


12 


52 


54 


5 


-29:27:17.1 


I 


188 


2 






— — 1.23 


79 


4 


92+* 


71 

58 




3C 280 


2210 


12 


56 


58 


2 


+47:20:22.2 


s 


45 


7 






— — 0.996 


44,8 


2 


25+* 


12 
96 




3C 294 


3207,3445 


14 


06 


44 





+ 34:11:26.1 


6,7 


190 


9 






— — 1.78 


76 


2 


90+* 


73 

60 


44.48 


WARP J1415. 1+3612 


4163 


14 


15 


11 





+36:12:03.5 


I 


88 


8 






— — 1.03 


80 


2 


82+* 


12 

96 




3C 324 


326 


15 


49 


48 


8 


+21:25:37.5 


S 


39 


4 






— — 1.21 


44 


4 


83+* 


28 
12 


43.82 


4C +15.55 


3229 


16 


25 


14 


4 


+ 15:45:22.8 


I 


51 









— — 1.406 


81 


3 


95+_l 


95 
78 


45.67 



Table 2: The final cluster sample, split by morphological class. 



Notes for Table [2] Nl - the NED gives z— 0.35, but the X-ray spectral data has a far better fit with z— 0.38. N2 - this cluster has two X-ray point sources 
within 25kpc, which overlap slightly. The second, at lOkpc, has log luminosity '^43.72. N3 - the NED gives two redshifts, but only one is given in the paper 
which NED refers to. Columns are: NED Name - Name of cluster in the NED. If more then one cluster name exists then the nearest is given. If there is 
no cluster within 2' then the nearest object name at the cluster redshift is given; Obs Ids - Chandra Observation ID; RA & DEC ~ Position of cluster as 
determined from the X-ray emission {J2000); Array - ACIS Detector or CCDs used; Exp - Good exposure time in ksec, after filtering (average over the selected 
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Good 

Target Name Observation Observation OBSID Array QSO 

RA DEC (ksce) 



HS0017+2116 


00 


20 


10 


80 


-f 21:32:51.00 


3063 


10.0 


ACIS 


s 


Y 


3C9 


00 


20 


25 


20 


+ 15:40:53.00 


1595 


17.5 


ACIS 


s 


Y 


WHDF 


00 


22 


33 


30 


+00:20:55.00 


2252 


71.0 


ACIS 


I 




GSGP4X:048 


00 


57 


17 


10 


-27:21:47.00 


2242 


10.5 


ACIS 


s 


Yt 


XMMlHR-3&:4 


01 


45 


38 


18 


-04:41:24.48 


4275,4276 


52.1 


ACIS 


I 




CADISOIHFIELD 


01 


47 


36 


20 


+02:20:03.30 


2240 


28.3 


ACIS 


I 




J0305+3525 


03 


05 


47 


40 


+35:25:13.40 


4142 


12.3 


ACIS 


s 


Y 


EXTENDEDCDF-S3 


03 


31 


48 


79 


-27:57:08.10 


5019,5020 


240.2 


ACIS 


I 




EXTENDEDCDF-S2 


03 


31 


52 


60 


-27:41:44.92 


5017,5018 


219.2 


ACIS 


I 




EXTENDEDCDF-S4 


03 


33 


01 


78 


-27:57:09.61 


5022 


78.7 


ACIS 


I 




EXTENDEDCDF-Sl 


03 


33 


06 


10 


-27:40:53.50 


5015,5016 


237.6 


ACIS 


I 




0406-244 


04 


08 


51 


50 


-24:18:16.50 


3058 


18.2 


ACIS 


s 


Y 


HS0818+1227 


08 


21 


39 


10 


+ 12:17:29.00 


3571 


19.7 


ACIS 


s 


Y 


0828+193 


08 


30 


53 


40 


+ 19:13:15.60 


3059 


17.4 


ACIS 


s 


Y 


APM08279+5255 


08 


31 


41 


60 


+52:45:16.80 


2979 


88.3 


ACIS 


s 


Y 


SDSS091316+ 591921 


09 


13 


16 


60 


+59:19:21.50 


3034 


9.8 


ACIS 


s 


Y 


QSO0910+564 


09 


14 


39 


30 


+56:13:21.00 


4821 


22.9 


ACIS 


s 


Y 


BRI0952-0115 


09 


55 


00 


10 


-01:30:05.00 


5194 


19.8 


ACIS 


s 


Y 


PClOOO+4751 


10 


03 


52 


80 


+47:36:54.30 


4152 


13.7 


ACIS 


s 


Y 


FSC10214+4724 


10 


24 


34 


50 


+47:09:09.80 


4807 


21.4 


ACIS 


s 


Y 


LH-NW-4 


10 


32 


06 


00 


+57:37:24.99 


3345 


38.3 


ACIS 


I 




LH-NW-6 


10 


33 


22 


00 


+57:55:25.00 


3343 


33.5 


ACIS 


I 




LH-NW-5 


10 


34 


02 


10 


+57:28:25.00 


3346 


38.1 


ACIS 


I 




LH-NW-9 


10 


35 


16 


00 


+57:46:24.99 


3348 


39.4 


ACIS 


I 




PC_1035+4747 


10 


38 


08 


20 


+47:31:36.60 


4154 


8.8 


ACIS 


s 


Y 


SWIRELOCKMAN7 


10 


43 


27 


23 


+59:10:15.07 


5029 


70.8 


ACIS 


I 




SWIRELOCKMANl 


10 


44 


46 


15 


+58:41:55.45 


5024 


63.7 


ACIS 


I 




SWIRELOCKMAN9 


10 


47 


13 


85 


+59:20:06.95 


5031 


65.0 


ACIS 


I 




SWIRELOCKMAN3 


10 


48 


32 


77 


+58:51:47.33 


5026 


68.7 


ACIS 


I 




Q1208+1011 


12 


10 


56 


90 


+09:54:26.80 


3570 


10.0 


ACIS 


s 


Y 


HDF-N 


12 


36 


49 


40 


+62:12:58.00 


2421,3293 


222.0 


ACIS 


I 




SDSSJ130216+003032 


13 


02 


16 


10 


+00:30:32.10 


3958 


10.7 


ACIS 


s 


Y 


SDSS1306+0356JE 


13 


06 


09 


30 


+03:56:43.50 


3966 


117.6 


ACIS 


s 


Y 


F864X:052 


13 


44 


07 


30 


-00:28:33.00 


2250 


9.5 


ACIS 


s 


Yt 


GROTH-WESTPHAL 


14 


17 


43 


60 


+52:28:41.20 


3305,4357,4365 


191.3 


ACIS 


I 




EGS-3 


14 


20 


28 


00 


+53:02:01.30 


5845,5846 


97.6 


ACIS 


I 




EGS-1 


14 


22 


42 


30 


+53:25:37.51 


5841,5842 


90.6 


ACIS 


I 




SDSSJ144231+011055 


14 


42 


31 


70 


+ 01:10:55.30 


3960 


10.8 


ACIS 


s 


Y 


DADDIFIELD 


14 


49 


09 


10 


+09:01:36.00 


5032,5033,5034 


87.2 


ACIS 


I 




QSO1508+5714 


15 


10 


02 


90 


+57:02:43.40 


2241 


88.5 


ACIS 


s 


Y 


ELAIS:N1 


16 


10 


21 


90 


+54:33:36.00 


888 


71.9 


ACIS 


I 




ELAIS:N2 


16 


36 


48 


48 


+41:01:45.90 


887 


73.1 


ACIS 


I 




2036-254 


20 


39 


24 


50 


-25:14:30.40 


3060 


19.6 


ACIS 


s 


Y 


2048-272 


20 


51 


03 


40 


-27:03:04.60 


3061 


17.7 


ACIS 


s 


Y 



Table 3. Observations of blank fields. The fields containing deliberately targeted QSOs are indicated - those marked f were in the 
"Extragalactic diffuse emission and surveys" category as a search for ROSAT identified NELGs, but no emission was seen at the target 
point so they were treated as blank fields. 



chips); Fx - Cluster approximate observed frame 0.5-8keV flux (10~ erg/cm /sec). All fluxes and luminosities are approximate as tire detector response 
was only calculated for the central pixel; L^l - Cluster approximate rest frame 0.5-8keV luminosity (10 erg/sec); Zj^2 - Cluster approximate rest frame 
0.1-2.4keV luminosity (10"^"^ erg/sec); rx'^ - reduced x'^ of the XSPEC fit to the cluster spectrum; z- Cluster redshift; Ref ~ Source of cluster redshift as listed 
below; Excess - Excess number of sources in the central 1 Mpc (compared to the prediction) with Icr errors; Centre - log Luminosity (in erg/sec) of any source 
detected within 25kpc of the cluster centre. 

Jll- lLaRoque et al.l |2003|); [ 21- IStocke et al.l )l99ll); f31- lEbeling et al.l i2007l); f 41- ISmail et"all il99^ ) ; [5]- iBohringer et ahl )2000l) : [61- | Vikhlinin et al.l 

I1998D : f71- IStruble fc RoodI |1999|); f81-lMullis et al.l J2003'); [91-'PerIman et aL H2002D; [10]- iBurenin et al.N2006l): fl ll-lde Grandi et al.Ul99g) ; [ 12]- Wright, 
Abies & All en (1983); [131- lAbazaiian et al.l )2003l); [14]- Allen et al. I2004b |); [151- ICaccianiga et a l. (2000); [16]- lEbeling et al.l il998l ); [17]- iGioia et al.l 

, al.l Jl99di); [211- Gomez, Hu ghes & Birkinshaw (2000); [221- Aller, Aller 
3|); [25p'Schindler et al. ( 200ll ); [261- Molthagen, Wendk er & Briel (1997); 
[301 - |Rome r et al . IgOOQi); [31l- | Gioia fc Luppind )199^: [321-IUlmer et al] 

^ _ 1) : [36l- |Ponman et al .1 )1994|) : [ 37]- [ Schindler et al.l il995l); [38 1- 

Holden et aH I2OO2I) ; [401- "Schade. Barrientos ' 



)1998|) : ri81- iTozzi et all ISo O al); [ 191- iDe Grandi fc Molendi (2002); [20]- Wei et 
& H ughes (1992); [ 231- lAbaz aiian ct al. 12005); [241- Roukema & Baitlik il99 
[27]- | EbeIing et al.l )1996|); [281 - Allen et al. (1992); [29]- Bo hri nger et a lj 120041) . 

200a); [331- iBorgani fc Guzzd j200ll); [34]- iFischer et al.l llOgsi ); [351- lEbeling et al.l )200^, . 

al.l I2OO2I) ; [401- Schade, Barrient os & Lopez-Cruz (1997) ; [41[- IBohringer et al.l i2004lT 



ICohen fc KneibI 1200^ 
]43]- | EdEe et al.lT200a ) 



[391- 
[44]- 



Spinrad et al.l )198I 



[491- Abell, Corwin fc Olowin (1989 ); [50]- 

- I Add man-. 

[59]- Caretta et al. 12002); [60]- Colless 12001); 

[65]- Postman, Geller & Huchra (1988) ^ 

'Schwope et al. (2000): [711- lPierre et al. | [19971); [721- jMoIinari et al 
[76]- Hewitt & Bur bidgg il99ll ); [771- IStanford et al, ' 



(206e): [471- 
[521- 



[451- D ^c et al. ( 19961); [461 - [B onamente et al 

_ [Couch et al. I lOol) ; [511 - [Pahle et al. 12002) 

Adelma n-McCarthv et al.' f2007f); [551- lAdelman-McCarthy et al^ f2003); [561 

Colless 12001); [61]- Delia Coca et al. 12000) 

[661- ,De Propris et al. (2007); [07]- Sand et al. (2005) 

11994): [731- Andreon et al. 12005); [741- 
Pentericci et al.iL2000.) ; [791- .Blakeslee et al 



[42]- [ 



Liang et al 



Giommi et i 



rj2000l); 



iHewitt fc Burbidgj jl989l) 



1200 



3);[78]- 



al.N2005|); [4 8]- ' Abazaiian et al.l 

Rizza et al.l il998l ); [531- Isarr icntos et al. ( 200-|); 
Jones et al. 12003); [57]- C ao, Wei fc H u (199 9) ; [58]- Henry et ^ 
[62]- White I200(j); [63]- iTucker et al.l JT991) ; [64]- Wittman et aQ 
[681-'Willick et al.l |2001D; [69]-|Arnaud et al. (.199^)^ 
Stern et al. |20()3); [75]- Delt orn et al.l 
12003.) ; [801- iElIis & Jonea >20(M|) ; [811- 



Cluster Name 
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RA (J2000) DEC (J2000) Net Counts Fx Sig HR (7„r 



3C. 


.295 


CXOGBA 


J141209.9- 


h520419 


14:12: 


:09, 


,9 


+52: 


:04: 


:19, 


,7 


70, 


,5 


36, 


,51 


16, 


,1 


-0.55 


0, 


,22 


3C. 


.295 


CXOGBA 


J141127.3- 


h521131 


14:11: 


:27, 


.4 


+52: 


;11: 


:31, 


,9 


73, 


,4 


21 


,42 


34, 


,0 


-0.19 


0, 


15 


3C. 


.295 


CXOGBA 


J141135.4- 


h521008 


14:11: 


:35, 


.4 


+52: 


:10: 


:08, 


,7 


10, 


,7 


4, 


,47 


5, 


,2 


1.00 


0, 


,00 


SC. 


.295 


CXOGBA 


J141132.9- 


h521103 


14:11: 


:33, 


.0 


+52: 


;11: 


:03, 


,5 


9, 


,7 


3, 


,95 


4, 


,9 


-1.00 


0, 


,00 


3C. 


.295 


CXOGBA 


J141132.2+521116 


14:11: 


:32, 


,3 


+52: 


:11: 


:16, 


,8 


7, 


,7 


2, 


,27 


3, 


,8 


-1.00 


0, 


,00 


3C. 


.295 


CXOGBA 


J141125.4+521047 


14:11: 


:25, 


,4 


+52: 


:10: 


:47, 


,4 


5, 


,8 


1, 


,70 


3, 


,0 


0.00 


1, 


,00 


SC. 


.295 


CXOGBA 


J141114.4- 


h520631 


14:11: 


:14, 


.4 


+52: 


;06: 


:31, 


,1 


60, 


,4 


26, 


,16 


21, 


,1 


-0.56 


0, 


,22 


3C. 


.295 


CXOGBA 


J141057.4- 


h521131 


14:10: 


:57, 


.4 


+52: 


:11: 


:31, 


,1 


27, 


,5 


8, 


,16 


13, 


,0 


-1.00 


0, 


,00 


3C. 


.295 


CXOGBA 


J141153.0- 


h521019 


14:11: 


:53, 


,1 


+52: 


:10: 


:19, 


,9 


25, 


,7 


11, 


,29 


10, 


,5 


-0.38 


0, 


,33 


SC. 


.295 


CXOGBA 


J141148.3- 


h521128 


14:11: 


:48, 


.4 


+52: 


;11: 


:28, 


,8 


10, 


,3 


4, 


,37 


4, 


,7 


-1.00 


0, 


,00 


3C. 


.295 


CXOGBA 


J141157.8- 


h520626 


14:11: 


:57, 


.8 


+52: 


;06: 


:26, 


,2 


30, 


,6 


14, 


,29 


9, 


,9 


-1.00 


0, 


,00 


3C. 


.295 


CXOGBA 


J141157.3- 


h520914 


14:11: 


:57, 


.4 


+52: 


:09: 


:14, 


,0 


9, 


,8 


4, 


,30 


3, 


,5 


0.00 


1, 


,00 


3C. 


.295 


CXOGBA 


J141123.4- 


h521332 


14:11: 


:23, 


,4 


+52: 


:13: 


:32, 


,1 


407. 


4 


118, 


,61 


188, 


,7 


-0.60 


0, 


,08 


3C. 


.295 


CXOGBA 


J141120.4+521210 


14:11: 


:20, 


,5 


+52: 


:12: 


:10, 


,3 


236. 


1 


68, 


,08 


60, 


,0 


-0.25 


0, 


,08 


3C. 


.295 


CXOGBA 


J141120.4- 


h521211 


14:11: 


:20, 


.4 


+52: 


;12: 


:11, 


,8 


208. 


8 


60, 


,19 


52, 


,3 


-0.75 


0, 


13 



Table 4. Properties of sources detected in the cluster fields. The full sample of sources in 150 cluster fields is given in the online journal 
and at http:/ /www. sc. eso.org/~rgilmour/. Columns are (1) cluster name given in Table [2l (2) full name of source (3)&(4) J2000 position, 
(5) net number of counts (6) flux in the 0.5-8keV band, xlO~^^erg/cm^/sec (assuming a spectrum with F = 1.7), (7) source significance 
as defined in Equation [7] (8) counts hardness ratio, (H-S)/(H+S), where H(2-8keV) and S(0.5-2keV) are set to for detections with 
significance < 3 in that band (9) error on hardness ratio, where good detections exist in both bands. 
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APPENDIX: TESTS FOR SYSTEMATIC 
ERRORS IN THE PIPELINE 

A.l Montecarlo simulations 

Small fluctuations in the background level, especially in the 
regions of intra-cluster emission, could cause sources to be 
missed by the wavelet detection method. In order to attempt 
to evaluate the detection efficiency of wavdetect near the 
flux limit, montecarlo simulations of faint sources were per- 
formed on cluster and blank field images. The difference in 
the number of sources detected, and with significance > 3, 
could then be evaluated as a function of radial position for 
the cluster and blank field samples. 

False sources were placed in 30 cluster and 30 blank 
field images, with up to 110 sources per image. The counts 
for each source corresponded to a multiple (1, 1.25, 1.5 etc.) 
of the fiux limit at the corresponding position on the image. 
This procedure is not straightforward, and may not produce 
accurate results as it is difficult to simulate X-ray sources, es- 
pecially those with very few photons, because of the complex 
nature of the Chandra PSF. The faint sources used in the 
montecarlo simulations were extracted from bright sources 
at the same off-axis radius, which accounts for the off-axis 
radial variation in PSF, but not any angular variation or 
difference between ACIS-I and ACIS-S detectors. The small 
number of suitable bright sources in the sample restricted 
the possible off-axis radii at which the false sources could 
lie. 

There is a small deviation between the detection rates 
in cluster and blank fields in the central 100 arcsec only. Sur- 
prisingly rather than the detection rate decreasing for the 
cluster fields, it increases slightly for the blank fields. This 
is most likely due to the problems with producing accurate 
input sources, as described above. Sources at the flux limit 
were around 15% more likely to be detected in the blank 
fields in the central 100 arcsec, but this difference decreases 
rapidly as the source flux increases. Combining the results 
for the cluster and blank fields with the Log N(> S) - Log 
S distribution gives an estimate of the number of sources 
missed in the cluster fields relative to the blank field predic- 
tion. This is found to be ~ 0.12 sources per cluster field in 
the central 25" — 100". Very few sources are expected to be 
missed in the central 25" due to the low area and high fiux 
limit. The maximum errors due to missed sources result in 
a ;^ lo" change in the results. 
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Figure 11. The average ratio of actual to expected (Poissonian) 
variance in N(> 5)for sources detected in the 0.5-8 keV band 
in 1000 sub samples of blank fields. Sub-samples of 1 - 5 fields 
arc investigated at three fiux values. This result has two limiting 
factors: Firstly the subsamples were picked from the same parent 
population (in particular at A'^ > 10~^^erg cm"'^ s~^ only 17 
blank fields were available) so the overlap between subsamples 
reduces the observed variance for moderate sample sizes. Secondly 
the values can be < 1 as the Poissonian errors were calculated 
from the mean field value, whereas in reality some fields are larger 
so have more sources and smaller errors. 



sources is small and the Poissonian errors completely explain 
the variance between the sub-samples. At 10~^*'^erg cm~^ 
s~^ the variation in individual fields is 20% larger than that 
expected from the Poissonian errors, which is attributable 
to AGN clustering. This effect decreases as the sample size 
is increased, and for samples of 5 fields the variation is only 
slightly above the expected value; samples of this size are 
therefore sufficient to largely average out the effects of large 
scale structure. At 10~^^erg cm~^ s"'^ there are few blank 
fields available so only the variance for individual fields is 
shown. This is ~ 10% higher than expected from Poissonian 
analysis, indicating that there is some effect due to large 
scale structure in faint sources too. In taking small samples 
of fields it is therefore advisable to use at least 5 fields in or- 
der to ensure that the stated errors are not underestimated 
due to clustering. 



A. 2 Cosmic variance in small samples 

The error calculations used here are based on the Poissonian 
errors on the sources detected, assuming that they are ran- 
domly distributed. However large scale structure may give 
rise to significantly larger errors in small samples. The clus- 
tering of X-ray sources appears to be stronger in low fiux 
sources than in high flux, and can give rise to sign i ficant 
field to field varia tion s (see for example Yang et al] |2003| . 
iMuUis et al. | |2004| and iBasilakos e^alTl2005^ . To check the 
magnitude of this effect in this survey, subsamples of 1 to 5 
blank fields were chosen at random. Figure [TT] shows the ra- 
tio between the observed variation in the Log N(>S) - Log S 
and the expected variation from Poissonian errors, for 1000 
subsamples of each size. At 10~^^erg cm~'^ s~^ the number of 



A. 3 Blank field results 

The results for the 44 blank fields were checked to ensure 
that the method and pipeline worked correctly. The blank 
field Log N(>S) - Log S was compared to the literature, and 
the radial distribution of blank field sources was compared 
to the pipeline prediction for the same fields. 

The blank field Log N(> S) - Log S dist r ibutio n was 
compared with that derived bv [Manners et all (|2003l l from 
the ELAIS fields, as shown in Figure [T2l and agrees to well 
within the la error bars. As with most blank field surveys, 
the ELAIS sample will be affected by sample variance as it 
only covers 2 Chandra fields. Unfortunately, all other blank 
field Chandra surveys calculate the Log N(> 5") - Log S for 
the 0.5-2 keV and 2-8 keV bands independently, so cannot 
be compared to the pipeline results directly. 
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Figure 12. The blank field Lo g N(>S) - Log S plot compared 
to that from the ELAIS fields of lManners et al.l | |2003|) . The solid 
line is the data from this paper. The crosses mark the data points 
and the dashed lines the Icr errors from the ELAIS fields. 



The combined radial distribution for all blank field 
sources is shown in Figure fTST al. In the galaxy cluster ob- 
servations the cluster is not generally placed at the centre 
of the detector array. Therefore, in order to reproduce the 
method used for the cluster fields, the radial distribution for 
the blank fields was measured from a point on the detector 
which corresponded to the cluster centre, in a randomly se- 
lected cluster observation. There is a slight excess of sources 
in the central 200" compared to the prediction. Although this 
is not very significant it does correspond to around one extra 
source per average field. However, the blank field sample is 
dominated by the four deepest observations, which actually 
only cover two regions of sky - the Chandra/Hubble Deep 
Field North and the Extended Chandra Deep Field South. 
These two regions are a factor of two deeper in exposure time 
than the average blank fields, and account for 13 per cent 
of the total sources, and 35 per cent of the fainter sources 
(< 10^^^'^erg cm^^sec^^). Because these two regions con- 
tain so many sources, a small variation in the number of 
sources carries more weight than for the shallower blank 
fields. The effect of cosmic variance is therefore amplified in 
Figure [T5r a), as explained in Section A. 2. 

Figure fTSf b) shows the same radial distribution, but ex- 
cluding the two regions of sky covered by the four deepest 
observations, to limit the effects of cosmic variance. The re- 
maining blank fields have far less range in exposure times, 
and as expected the excess in the radial distribution seen in 
Figure llSf a) disappears. This sample of blank fields has a 
similar range of exposure times to the cluster sample. Figure 
I13f b) therefore shows that the prediction is accurate so long 
as the source counts are not dominated by a few fields. Sec- 
tion A. 2 gives a more quantitative analysis of this problem. 

Figure fTST c) shows the radial distribution for all sources 
brighter than 10~^''erg cm~^sec~^. Again the source distri- 
bution matches the prediction well. At this fiux level all 
blank fields contribute to the source counts, and there is no 
problem due to cosmic variance. Figures flST b) and (c) can 
be compared to Figure lO^ a)- (c) to show that the excess seen 
in the cluster fields is real. 



2.0 



1.5 



-0.5 



-1.0 



2.0 




Sources ^/ 
Combined 1 -Sigma error i 



100 200 300 400 

Distance from centre (arcsec) 



500 



(a) 



1.5 - 



Sources 

Combined 1 -Sigma error 




100 200 300 400 

Distance from centre (arcsec) 



500 



(b) 



2.0 



1.5 



1.0 



0.5 



Sources 

Combined 1 -Sigma error 




100 200 300 400 

Distance from centre (arcsec) 



500 



(c) 



Figure 13. The radial prediction and actual distribution of 
sources in blank fields, with Icr errors, for (a) all sources in all 44 
blank fields; (b) all sources in 40 blank fields, excluding the four 
fields with very deep exposures, which only cover two regions of 
the sky. These two regions bias the sample as they contribute a 
significant fraction of the sources and are both overdense due to 
cosmic variance; (c) sources brighter than 10~^*erg cm^'^sec"^, 
in all 44 blank fields. 
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Figure 14. Fractional diflference between the Log N(>S) - Log 
S plots for ACIS-S and ACIS-I blank fields. The la errors are 
shown relative to the ACIS-I line, and are found by combining the 
errors on the two Log N(>S) - Log S distributions. The ACIS-S 
distribution is around Icr lower than the ACIS-I distribution at 
all fluxes. 



Figure 15. The difference between the total number of sources 
> 10~^''erg cm~^sec~^ when the 44 blank fields are split into 
two equal subsamples, for 1000 randomly chosen subsamples. The 
thick vertical line shows the value when the sample is split into 
the 22 ACIS-I and 22 ACIS-S fields (see Figure [Tit. This value 
is consistent with 22 randomly chosen fields (the probability of 
having two such different values randomly is 37 per cent) and 
so there is no evidence of a systematic offset in the number of 
sources detected in ACIS-I and ACIS-S blank fields. 



As a final check the total excess or deficit of sources 
(over the full radius) compared to the prediction, was ex- 
amined for each blank field. There was no noticable correla- 
tion between the deviation from the prediction and factors 
such as exposure time, or whether a field was merged or not. 
There was a small and insignificant correlation with ACIS 
array, as explained in Section A. 4 below. 



A. 4 The CCD array 

To check for systematic offsets between the pipeline results 
for fields observed with the ACIS-I and ACIS-S detectors, 
the radial distribution and Log N(> 5") - Log S distributions 
for blank fields observed with each detector were compared. 
It is also desirable to check for differences between the true 
blank fields and those which targeted high redshift QSOs. 
Unfortunately these cannot be done independently as the 22 
high redshift QSO fields were all observed with ACIS-S, and 
the 22 true blank fields were observed with ACIS-I. 

Figure [TD shows the difference between the Log N(>S) - 
Log S distributions for the 'true' blank and QSO fields. The 
QSO fields (ACIS-S) have a Log N(>S) - Log S distribution 
that is around lu lower than the 'true' blanks (observed 
with ACIS-I), so there are no significant extra sources in 
the high redshift QSO fields, and they are valid blank fields. 

The Icr offset between the ACIS-I and ACIS-S (which 
correspond to the QSO and 'true' blank fields) is of low sig- 
nificance, but it is worth checking that it is not a systematic 
error. The small difference in the source size between ACIS-I 
and ACIS-S images cannot account for the la variation (Sec- 
tion I5.1|l . As the ofi'set is also seen at higher fluxes, where 
the sky area is given by the total area of the detector, errors 
in the calculation of Smin also cannot explain the difference. 
In addition the flux calibration between ACIS-I and ACIS-S 



is accurate to within < 



whereas a 10% offset would be 



required to change the Log N(>S) - Log S by \a. Finally, 
the effect of sources overlapping at the edges of the images 
(ACIS-I observations have more large sources) is minimal, 
even for the deepest flelds. 

Instead, the small offset between ACIS-I and ACIS-S 
number counts is most likely to be just due to the high 
statistical variance between the flelds. The 44 blank fields 
were split into two equal subsamples and the difference in 
N{S > 10~^''erg cm~^ s~^) was computed. This was re- 
peated for 1000 randomly chosen subsamples. Figure 1151 
shows that the difference between the ACIS-I and ACIS- 
S samples is fully consistent with randomly chosen samples 
of blank fields. The chance of getting a difference of > 23 
sources between the two samples is ~37 per cent, which is 
in full agreement with the size of the la error bars in Figure 

M 

The radial distributions for the ACIS-I and ACIS-S 
fields were also compared. Comparing the results for ACIS- 
I and ACIS-S blank fields against the prediction from the 
blank field Log N(>S) - Log S, the ACIS-S flelds end la be- 
low the prediction in agreement with Figure [HI Otherwise 
both distributions are fiat to within the errors, and there is 
no link between over or under prediction and radius. 



^ see details in 

http: / /cxc. harvard.edu/cal/docs/cal_present_status. html 



